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ABSTRACT 


This  report's  basic  task,  which  has  been  successfully 
accomplished,  is  that  of  developing  the  mathematical  model 
simulating  the  following  cargo  unloading  procedure. 

The  cargo  is  to  be  unloaded  from  a  ship  at  some  distance 
from  the  shore.  The  final  destination  of  the  cargo  is  to  be 
at  point  A  on  the  beach  where  no  port  facilities  are  available. 
The  transferring  of  cargo  from  the  ship  to  the  shore  is  to  be 
accomplished  by  means  of  transfer  vehicles,  such  as  amphibious 
craft.  The  cargo  transfer  from  the  ship  to  the  transfer 
vehicles  is  to  be  accomplished  by  ship-based  unloading  gear, 
such  as  ship-based  cranes.  The  cargo  transfer  from  the  trans¬ 
fer  vehicles  on  the  shore  is  to  be  accomplished  by  beach-based 
unloading  gear,  such  as  fork  lifts. 

The  above  mentioned  simulation  enables  the  user  to  gain 
insight  into  the  aforementioned  unloading  procedure  and  thus 
to  derive  correctly  the  optimal  use  strategy  of  the  ship-to- 
shore  transfer  operations,  which  is  the  ultimate  goal  of  the 
present  analysis. 
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1 .  Introduction 


The  basic  task  of  the  problem  posed  ic  that  of  develop¬ 
ing  the  methodology  that  will  permit  the  overall  analysis  of 
a  cargo  offloading  procedure.  In  the  offloading  procedure 
under  investigation,  the  cargo  is  to  be  unloaded  from  a  ship, 
henceforth  referred  to  as  the  mother  ship,  which  is  at  some 
distance,  say  x  miles,  from  the  shore.  The  Mother  Ship  is  to 
be  stationary  during  the  entire  unloading  operation.  The 
final  and  only  destination  of  the  cargo  is  to  be  a  point  A  on 
the  beach.  There  are  to  be  no  port  facilities  on  the  shore 
or  beach  unloading  areas. 

The  cargo  involved  in  this  study  is  to  be  contained  in 

i)  Containers  or  pallets  of  arbitrary  size,  weight 
and  capacity  that  are  not  capable  of  any  self- 
induced  motion,  or 

ii)  Vehicles  of  arbitrary  size,  weight  and  capacity 
that  are  capable  of  self-induced  rolling  motion 
only.  It  should  be  noted  that  in  this  case  the 
vehicle  itself  may  be  the  cargo. 

The  actual  transferring  from  the  Mother  Ship  to  the 
beach  is  to  be  accomplished  by  means  of  amphibious  craft, 
whose  number  and  characteristics  hu  <e  been  prespecified. 

These  are  henceforth  referred  to  as  the  transfer  vehicles, 
such  as  LARCs,  GEMs,  etc.  The  loading  into  the  transfer 
vehicles  alongside  the  mother  ship  is  to  be  ir  one  of  the  two 


modes:  sequential  or  simultaneous.  Similarly,  but  toally 
independent  from  the  loading  mode,  the  unloading  from  the 
transfer  vehicles  on  the  beach  is  to  be  in  one  of  the  two 
modes:  sequential  or  simultaneous . 

* 

In  order  to  accomplisn  the  cargo  transfer  from  the  Mother 
Ship  to  the  transfer  vehicles,  the  Mother  Ship  is  to  be  pro¬ 
vided  with  all  the  necessary  unloading  facilities,  for  example 
ship-borne  crane (s),  ramp(s),  etc.  However,  the  transfer 
vehicles  are  not  to  be  provided  with  any  special  unloading 
facilities,  because  some  means,  such  as  a  fork  lift,  is  to  be 
made  available  on  the  beach  tc  carry  out  the  cargo  unloading. 

Finally,  the  prespecified  number  of  transfer  vehicles 
and  beach  unloading  gear  is  to  be  made  available  at  point  A, 
and  their  transportation  and  arrival  is  to  be  independent  of 
that  of  the  Mother  Ship,  and  from  each  other. 

With  the  above  description,  the  cargo  offloading 
procedure  under  investigation  has  been  fully  defined.  In 
order  to  complete  the  description  of  the  problem  posed,  it 
remains  to  define  the  analysis  objectives,  which  can  be 
stated  as  follows:  The  resulting  technique  is  to  be  designed 
to  first  provide  a  common  measure  of  success  for  a  number  of 
prespecified  use  strategies  for  given  ship-based  loading 
facilities,  transfer  vehicles  and  beach-based  unloading 
facilities  distributions,  for  a  given  x  and  environment  state, 
and  for  given  breakdown  considerations.  The  common  measure 
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of  success  is  to  depend  on  time  and/or  level  of  risk  or  un¬ 
certainty.  Thus  the  final  analysis  objective  is  to  determine 
the  optimum  strategy  (minimum  time  and/or  level  of  risk  or 
uncertainty)  among  those  examined  or,  if  the  findings  of  the 
previous  calculations  suggest  it,  to  continue  the  analysis 
with  new  strategies  until  the  optimal  one  is  found. 

With  the  above,  the  problem  description  has  been  com¬ 
pleted.  The  solution  process  is  to  be  outlined  in  the  sub¬ 
sequent  Sections  and  Appendices.  We  start  our  discussion 
with  a  general  outlii  e  of  the  solution  process. 


2.  Solution  Process 


The  solution  process  in  this  study  is  best  illustrated  by 
the  block  diagram  shown  in  Fig.  2-1,  A  brief  discussion  of 
each  step  involved  in  the  solution  process  is  given  below. 

Problem  Definition.  This  involved  the  following: 

i)  Identification  of  this  variables  of  the  problem. 

ii)  Establishment  of  the  relations  amcng  the 
variables. 

iii)  Identification  of  the  dependent  and  independent 
variables  of  the  problem  consistent  with  i) 
and  ii ) . 

iv)  Establishment  of  the  range  of  variation  of  all 
problem  parameters. 

v)  Selection  of  the  figure  of  merit. 

By  definition,  the  variables  of  a  problem  are  those 
parameters  necessary  to  fully  describe  a  given  system  to  the 
degree  of  accuracy  and  extent  desired.  The  process  of  var¬ 
iable  identification  for  a  new  problem,  such  as  ours,  is  a 
major  and  a  very  difficult  task.  In  order  to  simplify  our 
task  of  identifying  our  variables,  it  was  found  advantageous 
to  first  identify  the  subsystems  involved  in  our  study  and 
then  find  the  variables  necessary  to  fully  describe  each  sub¬ 
system  to  the  degree  of  accuracy  and  extent  desired.  The  sub¬ 
systems  involved  in  our  study  were  found  to  be: 
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1.  The  mother  ship. 

2.  The  payload. 

3.  The  ship-based  unloading  facilities,  S.U.F. 

4.  The  shore-based  unloading  facilities,  B.U.F. 

5.  The  transfer  vehicles,  T.V. 

The  variable  identification  is  deferred  until  the  next 
Section  where  the  Problem  Definition  will  be  discussed  in  de¬ 
tail.  It  is  of  importance  to  note  in  conjunction  with  the 
discussion  of  the  variable  identification  that  although  it  was 
recognized  that  the  environment  state  as  described  by 

1.  Wind  speed, 

2.  Sea  state, 

3.  Current, 

4.  Tide, 

5.  Obstacles, 

6.  Beach  configuration,  and 

7.  Shore  configuration 

influences  our  offloading  operation,  it  was  decided  to  de¬ 
scribe  their  effect  by  externally  adjusting  the  magnitude  of 
the  appropriate  parameters.  For  this  reason,  it  is  not  neces¬ 
sary  to  identify  any  variables  that  describe  the  environment 
state.  However,  when  assigning  the  i^agnitude  of  the 
parameters  that  are  affected  by  the  environment  state,-  the 
user  must  assign  the  appropriate  values  by  taking  the  environ¬ 
ment  state  into  consideration. 


The  establishment  of  the  relations,  if  any,  among  the 
problem  variables  is  a  necessity  as  they  form  a  part  of  the 
mathematical  model.  The  reason  for  this  is  that  because  the 
mathematical  model  is,  by  definition,  the  replica  of  our  sys¬ 
tem  in  the  form  of  mathematical  equations,  any  such  relations 
must  be  part  of  it  in  order  to  permit  it  to  be  a  true  replica 
of  the  original.  The  establishment  of  any  such  relations  is 
deferred  until  the  next  Section  where  the  Problem  Definition 
will  be  discussed  in  detail. 

As  was  stated  at  the  beginning  of  this  section,  the 
problem  variables  were  so  selected  that  when  the  appropriate 
values  were  assigned  to  them,  they  could  serve  to  define  a 
given  offloading  system  uniquely  to  the  degree  of  accuracy  and 
extent  desired.  This,  however,  should  not  be  taken  to  imply 
that  if  we  assign  arbitrary  but  logical  values  to  these  var¬ 
iables  we  will  always  be  able  to  generate  an  offloading  system 
because  of  the  possible  interrelations  among  the  variables, 
which  do  not  permit  independent  selection  of  values  for  all 
the  interrelated  variables.  This  fact  makes  it  necessary  to 
identify  those  variables  whose  values  can  be  assigned  arbitrar¬ 
ily  and  those  whose  value  cannot.  This  is  because  our  solu¬ 
tion  method  involves  the  evaluation  of  different  offloading 
systems  which  are  generated  by  prespecifying  the  values  of 
their  variables.  This  is  best  done  by  classifying  the  var¬ 
iables  into  dependent  and  independent  ones.  The  independent 
variables  are  those  variables  which  must  be  prescribed  to 


completely  describe  our  system  as  desired.  The  dependent 
variables  are  the  ones  remaining  in  our  original  list  of 
variables  after  the  independent  ones  have  been  removed.  The 
identification  of  the  dependent  and  independent  variables  is 
deferred  until  the  next  Section  where  the  Problem  Definition 
will  be  discussed  in  detail. 

Ideally  speaking,  we  would  like  to  impose  no  restric¬ 
tions  on  the  range  of  variation  of  our  problem  parameters,  as 
this  will  tend  to  decrease  the  universality  of  our  methodology, 
something  a  true  engineer  does  not  like  to  do  voluntarily. 
However,  there  are  practical  considerations,  common  to  this 
type  of  problem,  which  make  it  necessary  that  we  impose  limits 
on  the  range  of  variation  of  our  problem  parameters.  The  most 
common  of  these  considerations  (requiring  us  to  compromise  by 
imposing  limits  on  the  range  of  variation  of  our  problem 
parameters)  is  that  it  is  impossible  to  construct  the  math¬ 
ematical  model  valid  over  the  entire  range  of  variation  of  the 
problem  parameters.  In  the  few  times  that  it  is  possible  to 
construct  such  a  model,  it  again  becomes  necessary  to  restrict 
the  range  of  variation  to  simplify  the  model  and  make  it  a 
useful  engineering  tool.  Therefore,  the  above  indicates  the 
need  for  the  introduction  of  restrictions  as  unavoidable . 

These  restrictions,  of  course,  ought  to  be  introduced  with 
great  care.  Care  should  be  taken  because  we  do  net  wish  to 
reduce  the  universality  of  our  methodology  unnecessarily,  as 


we  wish  to  utilize  it  to  solve  most,  if  not  all,  of  the 
problems  that  we  are  likely  to  encounter,  but  at  the  same 
time  we  wish  to  obtain  this  solution  with  relative  ease  and 
consistent  but  adequate  precision.  The  introduction  of  these 
restrictions  is  deferred  until  the  next  Section,  when  the 
Problem  Definition  will  be  discussed  in  detail. 

Finally,  as  was  already  mentioned  in  the  Introduction, 
the  figure  of  merit  (the  measure  of  success)  is  the  weighted 
combination  of  time  and  level  of  risk  and  uncertainty.  The 
time  component  of  the  figure  of  merit  involves  the  calculation 
of  time  elapsed  since  the  start  of  the  mission  to 

j .  prepare  the  mother  ship  for  departure  after  all 
the  payload  has  been  transferred  into  the  transfer 
vehicles,  and  all  transfer  vehicles  have  cleared 
the  mother  ship, 

2.  complete  the  payload  transfer  from  the  mother  ship 
to  point  A  on  the  beach, 

3.  return  all  transfer  vehicles  to  their  bases,  and 

4.  return  all  beach-based  unloading  facilities  to 
their  bases. 

The  level  of  risk  and  uncertainty  component  of  the  figure  of 
merit  involves  the  calculation  of  the  percentage  of  transfer 
vehicles  that  did  not  complete  their  mission  because  of 


breakdown.  The  factors  determining  the  likelihood  of  break¬ 
down  for  each  transfer  vehicle  are 

1.  Reliability  of  the  transfer  vehicle's  components. 

2.  Hazard  vulnerability. 

3.  Control  stability. 

4.  Operational  limitations. 

The  establishment  of  the  exact  nature  of  the  figure  of  merit 
is  deferred  until  the  next  Section  where  the  Problem  Defini¬ 
tion  will  be  discussed  in  detail. 

Formulation  of  the  Mathematical  Model. 

As  was  stated  earlier,  the  mathematical  model  is,  by 
definition,  the  replica  of  our  system,  to  the  degree  of  accur¬ 
acy  and  extent  desired,  in  the  form  of  mathematical  equations. 
Because  of  the  nature  of  our  problem,  the  mathematical  model 
is  stochastic  in  nature.  In  setting  up  the  mathematical  model, 
care  was  taken  to  keep  it  as  simple  as  possible  to  permit  easy 
analysis,  and  yet  to  construct  it  so  that  it  exhibits  all  the 
phenomena  under  consideration,  as  required.  The  actual  con¬ 
struction  of  the  mathematical  model  is  deferred  until  the 
fourth  section. 

Selection  of  Solution  Method. 

From  the  Problem  Definition  one  may  easily  observe  that 
the  easiest  way  to  achieve  the  desired  goal,  namely,  to  find 
the  optimum  use  strategy,  is  to  treat  the  use  strategy  as  a 
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problem  variable,  and  then  solve  the  problem  under  investiga¬ 
tion  as  an  optimization  problem.  The  resulting  optimization 
problem  is  a  mixed  integer  one,  or  simply  an  integer  problem 
if  the  waiting  times  of  the  transfer  vehicles  and  unloading 
facilities  are  approximated  by  integer s.  However,  it  was  soon 
discovered  that  in  order  to  solve  the  problem  as  an  optimiza¬ 
tion  one,  drastic  simplifications  had  to  be  made  to  the  math¬ 
ematical  model  to  allow  us  to  efficiently  implement  the  solu¬ 
tion  in  present-day  computers.  The  reason  for  this  is  that 
all  the  solution  methods  required  most,  if  not  all,  of  the 
past  history  of  the  system  to  be  stored  in  the  computer 
memory.  The  drastic  simplifications  necessary  made  our  meth¬ 
odology  a  very  inefficient  engineering  tool.  Last,  but  not 
least,  if  the  use  strategy  obtained  by  solving  the  problem  as 
an  optimization  one  was  very  complicated,  it  probably  would 
have  been  very  difficult  to  implement  in  practice  (because  the 
system  might  be  operating  under  external  pressure) .  For  this 
reason  it  would  be  very  easy  to  vio. ate  a  complicated  optimum 
use  strategy  and  to  actually  adopt  a  suboptimal  solution 
whose  merit  cannot  be  estimated  in  any  way,  which  is  a  very 
undesirable  situation  as  it  defeats  the  purpose  of  this 
analysis. 

For  the  reasons  given  above  it  was  decided  to  develop  a 
methodology  that  will  yield  the  desired  solution,  not  neces¬ 
sarily  as  directly  as  it  would  have  been  provided  by  the 
optimization  theory,  but  one  which 


1.  could  be  implemented  without  requiring  major 
simplifications  in  the  mathematical  model,  and 

2.  could  test  logical  and  likely  optimal  use 
strategies  which  have  a  very  high  probability  of 
being  implemented  in  practice. 

The  method  that  satisfied  all  the  above  requirements  was 
the  digital  simulation  method,  which  was  utilized  in  this 
study  to  obtain  the  desired  solution.  By  this  method  the 
optimal  solution  was  obtained  by  testing  different  use  strat¬ 
egies  that  satisfied  the  second  requirement  given  above.  To  . 
develop  these  use  strategies,  one  is  guided  by  logic,  es¬ 
pecially  in  developing  tne  first  use  strategy  to  be  tested, 
and  by  the  insight  gained  from  the  previous  tries  when  this 
is  available.  In  theory  the  true  optimum  can  be  found  by  ex¬ 
amining  all  the  possible  use  strategies,  and  although  this 
can  be  a  large  number,  it  is  always  finite.  However,  from 
the  above  it  should  be  clear  that  this  is  not  necessary,  as 
we  wish  to  find  the  optimum  solution  that  has  also  a  very 
high  probability  of  being  implemented  in  practice.  Fortunate¬ 
ly,  this  can  oe  achieved  by  a  reasonably  small  number  of 
tries.  The  discussion  of  the  actual  details  of  the  method¬ 
ology  used  in  this  study  is  deferred  until  the  fifth  section. 


Problem  Solution 


This  involves  the  preparation  of  the  input  required  by 
tne  computer  program.  Special  care  must  be  taken  when  pre¬ 
paring  the  input  of  the  parameters,  whose  magnitude  is 
affected  by  environment  state  and  breakdown  considerations. 
Further  discussion  on  this  topic  is  deferred  until  later. 

Evaluation  of  Results 

With  reference  to  Fig.  2-1,  special  care  was  taken  that 
the  only  iteration  required  in  the  Solution  Process  is  the 
preparation  of  new  input  data  for  the  examination  of  a  new 
use  strategy,  if  the  Evaluation  of  Results  suggests  it.  It 
is  anticipated  that  there  will  never  be  any  need  for  the  al¬ 
teration  of  the  first  three  steps  of  Fig.  2-1,  as  care  was 
taken  to  make  the  methodology  developed  here  general,  in 
order  to  handle  all  cases  likely  to  be  encountered  in  prac¬ 
tice.  However,  if  a  case  arises  where  a  change  must  be  intro¬ 
duced  in  these  three  steps,  the  user  must  read  very  carefully 
the  next  three  sections  so  that  he  may  correctly  alter  the 
present  method  to  suit  nis  needs.  Further  discussion  on  this 
topic  is  deferred  until  later. 

The  above  completes  the  introduction  in  the  Solution 
Process >  In  the  next  section  a  detailed  presentation  of  the 
Problem  Definition  will  be  given. 


3.  Problem  Definition 

In  this  section  a  detailed  discussion  of  each  step  in¬ 
volved  in  the  Problem  Definition  is  given. 

i)  Identification  of  the  Problem  Variables 

In  view  of  the  fact  that  the  smaller  the  number  of 
variables  in  a  problem  the  more  economical,  and  in  many- in¬ 
stances  the  more  efficient,  the  solution  process  becomes,  an 
attempt  was  made  to  keep  the  number  of  variables  of  this 
problem  to  a  minimum.  To  do  so  it  was  necessary  to  introduce 
certain  assumptions.  However,  special  care  was  taken  so  that 
the  nature  and  number  of  these  assumptions  was  such  as  not  to 
diminish  the  generality  of  our  methodology .  These  assumptions 
will  be  enumerated  in  the  fourth  topic  of  this  section,  when 
the  range  of  variation  of  the  problem  parameters  is  discussed. 

As  was  mentioned  earlier,  in  order  to  simplify  our 
task  of  identifying  our  problem  variables,  the  subsystems 
involved  in  our  study  were  identified.  These  are  shown 
diagrammatical ly  in  Fig.  3-1.  For  presentation  purposes,  the 
variables  that  are  utilized  to  define  each  subsystem  and  at 
the  same  time  appear  in  the  computer  input  will  be  listed 
first,  while  the  remaining  variables  necessary  to  complete 
each  subsystem's  description  will  be  given  later. 
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Fia.  3-1  GENERAL  MODEL  CONSTRUCTION 


MOTHER  SHIP  DESCRIPTION 


The  variables*  selected  to  describe  the  Mother 

Ship's  performance  are: 

TAM  giving  the  number  of  units  of  time  after  the  start  of 
the  mission  that  the  Mother  Ship  is  expected  to  arrive 
in  the  theater  of  operations  [-999<TAM<9999  ('treated 
as  a  floating  point  number j]. 

T1  giving  the  number  of  units  of  time  required  to  com¬ 
plete  the  mooring  operations  of  the  Mother  Ship  after 
it  arrived  in  the  theater  of  operations  [0 . <T1<9999 
(treated  as  a  floating  point  number^]. 

T2  giving  the  number  of  units  of  time  required  to  free 
the  Mother  Ship  from  its  moorings  and  make  ready  to 
travel  after  all  S.U.F.  are  secured  to  position  and 
all  T.V.  have  cleared  the  Mother  Ship  [0 . <T2<9999 
(treated  as  a  floating  point  number^]. 

IDMA**  indicating  the  nature  of  the  process  concerning  the 
arrival  of  the  Mother  Shijx 

ID1 *  *  indicating  the  nature  of  the  process  concerning  the 
mooring  operation  of  the  Mother  Ship. 

ID2**  indicating  the  nature  of  the  process  concerning  the 

operation  of  freeing  the  Mother  Ship  from  its  moorings. 


•The  magnitude  of  all  the  variables  selected  to  describe 
the  Mother  Ship's  per;  jrmance  except  the  magnitude  of 
INMS,  INI ,  IN 2  and  TIME  is  affected  by  environmental  state. 

**IDMA,  ID1,  ID2  =2  if  the  process  is  deterministic. 

=  2  if  the  process  is  stochastic,  drawn 
from  U(0,1 )  distribution. 

=  3,4. . . 9999  if  the  process  is  stochastic, 
drawn  from  distributions  to  be  devel¬ 
oped  by  the  user,  if  so  desired. 


INMS  giving  the  seed  of  the  U(0,1 )  distribution,  which 
serves  to  predict  the  stochastic  behavior  of  the 
Mother  Ship's  arrival  process,  if  IDMA  =  2*** 
[1<INMS<_10*-1  ]  • 

AMINM  giving  the  minimum  value  of  the  range  of  variation  of 
TAM,  if  IDMA  =  2***1-999.99999  <  AMINM  <  9999.99999 ]. 

AMAXM  giving  the  maximum  value  of  the  range  of  variation  of 
TAM,  if  IDMA  -  2***1-999.99999  <  AMAXM  <  9999. 99999]. 

INI  giving  the  seed  of  the  U(0,1)  distribution,  which 

serves  to  predict  the  stochastic  behavior  of  the  Mother 
Ship's  mooring  operation,  if  ID1  =  2*  [1<IN1<109-1] . 

AMIN1  giving  the  minimum  value  of  the  range  of  variation  of 
T1  if  ID1  =  2* [0 . <AMIN1<9999 .  99999] 

AMAX1  giving  the  maximum  value  of  the  range  of  variation  of 
T1  if  ID1  =  2*[0. <AMAX1<9999.  99999] . 

IN2  giving  the  seed  of  the  U(0,1)  distribution,  which 

serves  to  predict  the  stochastic  behavior  of  the  opera 
tion  of  freeing  the  Mother  Ship  from  its  moorings, 
if  ID2  -  2**[1<ID2<10*-1) . 

AMIN2  giving  the  minimum  value  of  the  range  of  variation  of 
T2  if  ID2  =  2**[0. <AMIN2<9999. 99999] . 

AMAX2  giving  the  maximum  value  of  the  range  of  variation  of 
T2  if  ID2  -  2**[0.  <AMAX2<_9999. 99999]  .  ! 

’l-IME  giving  the  units  of  time  utilized  in  this  study 
[0  <  TIME  <  12  alpha  numeric  characters] . 


*If  ID1-1  INI,  AMIN1  and  AMAX1  are  not  problem  variables. 

**If  ID2*1  IN2,  AMIN2  and  AMAX2  are  not  problem  variables. 

*  * *If  IDMA-1  INMS,  AMINM  and  AMAXM  are  not  problem  variables. 


PAYLOAD  DESCRIPTION 


N 


The  variables  selected  to  describe  the  payload  are: 


k  giving  the  number  of  payload  units  to  be  unloaded  by  each 

of  the  k  (k  =  1,2. ..K)  S.U.F. 

K 

1  <  N.  <  1000  and  7  N.  <  1000 
*  k=l  K 


WCn  giving  the  weight  of  each  of  the  n 


K  ) 


n  *=  1,2.  .  .  £  N, 


k-1 


payload  units  plus  that  of  their  lashings  [0. <WCn<99.  99] . 


VCn  giving  the  volume  of  each  of  the  n 


n  «  1,2. . .  lM. 


fc-I 


payload  \mits  together  with  that  of  their  lashings* 
[0. <VCn<9999  (treated  ae  a  floating  point  number)]. 

WGHT  giving  the  units  of  weight  utilized  in  this  study 
[0<WGHT<12  alpha  numeric  characters] . 


VOL  giving  the  units  of  volume  utilized  in  this  study 
[0<V0L<12  alpha  numeric  characters] . 


*If  all  the  T.V.  employed  in  this  study  do  not  permit  vertical 

stowing  of  the  payload,  VC can  be  utilized  to  give  the  volume 

per  unit  height  plus  the  surface  area  required  by  their  lash- 

K 

ings  rather  than  the  volume  of  each  of  the  n  (n*l,2. . . 
payload  units  and  that  of  their  lashings,  as  this  is  a  quan¬ 
tity  much  easier  to  estimate. 


S.U.F.  DESCRIPTION 


The  variables*  selected  tc  describe  the  S.U.F.'s  per¬ 
formance  are: 


K  giving  the  number  of  S.U.F.**  involved  in  this 

case  [1  £  K  £  20] . 

ISUD***  indicating  the  nature  of  the  unloading  mode  at  the 
Mother  Ship. 

ISCSL****  indicating  which  of  the  S.U.F.  use  strategies  is  to 
be  used. 

TSCj  £  giving  the  time  required  for  the  fcth  (k  *=  1,2...  K) 
S.U.F.  to  be  made  ready  to  start  the  unloading  oper 
ation  and  reach  the  fcth  ship  unloading  area  after 
the  Mother  Ship  is  properly  moored  [ 0 .  <  TSC 1  ^ 

<  9999  (treated  as  a  floating  point  number)]. 


*  The  magnitude  of  all  the  variables  selected  to  describe  the 
S.U.F.'s  performance  except  the  magnitude  of  INSCa  ^ 

(j  -  1,2  ...  5,7, 8,9;  k  -  1,2  ...  K)  is  affected  2*by 
th|  environment  state. 

**  Each  S.U.F.  is  identified  by  a  distinct  number,  k,  such 
that  1  <  k  <_  K. 

***  ISUD  *  1  if  the  unloading  mode  at  the  Mother  Ship  is  to 
be  in  parallel. 

ISUD  *  2  if  the  unloading  mode  at  the  Mother  Ship  is  to 
be  sequential. 

****  jf  isCSL  ■  1,  S.U.F.  use  strategy  SLSCA  is  used  to  se¬ 
lect  the  appropriate  S.U.F.  when  necessary. 

•2,  S.U.F.  use  strategy  SLSCB  is  used  to  se¬ 
lect  the  appropriate  S.U.F.  when  necessary. 

* 3,4  ...  9,  additional  S.U.F.  use  strategies 
to  be  developed  by  the  user,  if  desired, 
for  selecting  the  appropriate  S.U.F.  when 
necessary. 
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y 

TSC0  .  giving  the  time  required  for  the  kth  (k  =  1,2... K) 

S.U.F.  to  travel  to  any  of  the  ^payload  units  from  the 
kth  ship  unloading  area  [ 0 .  <  TSC  ^  £  9999  (treated 
as  a  floating  point  number)]. 

TSC s  £  giving  the  time  required  to  release  any  of  the 

(k  =  1,2... K)  payload  units  after  the  kth  S.U.F. 
has  reached  the  payload  unit  in  question  [0.  <  TSC s  ^ 

<_  9999  (treated  as  a  floating  point  number)]  . 

TSCa  giving  the  time  required  to  secure  any  of  the 

(k  =  1,2... K)  payload  units  on  the  kth  S.U.F.  after 
the  payload  unit  in  question  has  been  released 
[0 .  <  TSC^  £  <_  9999  (treated  as  a  floating  point 
number) ] . 

TSC 5  £  giving  the  time  required  to  transport  any  of  the 

(k  -  1, 2. . . K)  payload  units  to  the  kth  ship  unloading 
area  after  the  payload  unit  in  question  has  been  se¬ 
cured  on  the  S.U.F.  [0 .  <  TSC „  .  <  9999  (treated  as  a 

o  j  K  — 

floating  point  number)] . 

TSC?  giving  the  time  required  to  unload  and  then  free  any 
of  the  (k  =  1,2...K)  payJoad  units  from  the  kth 
S.U.F.  and  to  then  make  the  kth  S.U.F.  ready  to 
travel  again.  This  operation  is  performed  only  if 
(a)  the  appropriate  T.V.  is  properly  secured  in  the 
kth  ship  unloading  area  and  has  completed  its  refuel¬ 
ing  (if  refueling  was  necessary) ,  (b)  the  previous 
payload  unit  unloaded  by  the  S.U.F.  in  question  la 
fully  secured  in  the  T.V.  in  question  (this  require¬ 
ment  is  void  if  the  payload  unit  in  question  is  the 
first  payload  unit  to  be  unloaded  in  any  of  the  T.V.'s 
trips)/  and  (c)  the  T.V. 's  remaining  capacity  can 
accept  the  payload  unit  in  question.  If  any  of  the 
above  is  not  satisfied,  the  kth  S.U.F.  must  wait 


TSC y  k  until  all  three  requirements  are  satisfied 

(cont'd  )  <  T^C ?  k  —  9999  (treated  as  a  floating  point 

number) ] . 

TSC &  k  giving  the  time  required  for  the  fcth  (k  =  1,2...K ) 
S.U.F.  to  travel  back  to  its  original  position 
from  the  kth  ship  unloading  area  after  the  last  of 
the  Nk  payload  units  has  been  transferred  onto  the 
appropriate  T.V.  [0.  <  TSC g  k  <  it 999  (treated  as  a 
floating  point  number)]. 

TSC 9  k  giving  the  time  required  for  the  kth  (k  =  1,2... K) 

S.U.F.  to  be  secured  to  its  original  position. 

IDSC  .  .*  indicating  the  nature  of  each  of  the  (3  =  1,2... 

3  3  *  ‘  ^ 

2  ...5,7, 8, 9)  processes  described  above  for  each  of 

the  k  (k  =  1, 2...K)  S.  U.F. 

INSC.  ,  giving  the  seed  of  the  U(0,1)  distribution,  which 
3%  i  * 

serves  to  predict  the  stochastic  behavior  of  each  of 
the  j  2  (q  2  -  1,2.  ..  5,7,8,  9)  processes  described  above 
for  each  k  (k  =  1,  2. .. K) ,  if  IDSC  ^  k  =  2** 

[1  <  INSC.  .  <  io 9  -1] . 

V2  *  * 

AMINSC .  k  giving  the  minimum  value  of  the  range  of  variation  of 

*  TSC k  (3Z  -  1,2.  ..  5,7 ,8,9 ;  k  «  1,2...K),  if 

IDSC  /  hm2**  [0.  <  AMINSC.  .  <  9999 .  99999]  . 

3  2  .»  *  J2  .»  * 

i4Af>l/SC  .  ^  giving  the  maximum  value  of  the  range  of  variation  of 
**  TSCj  k  (32  =  1,2.  ..  5,7,8,9s  k  =1,2...K),  if 
IPSc/  .  *  2**  [0.  <  AMAXSC  .  .  <  9999 .  99999]  .  ' 

( _ _  -  -  3*  .k -  3  it  k 

*IDSC  .  1,-1  if  the  process  is  deterministic. 

=2  if  the  process  is  stochastic,  drawn  from  U(0,1) 
distribution. 

=  3,4. . .9999  if  the  process  is  stochastic,  drawn  from 
distributions  to  be  developed  by  the  user,  if  so 
desired. 

**If  IDSC.  .*2  INSC.  .,  AMINSC .  .  and  AMAXSC.  .  are  not 

3 2>k  32jk*  32  ,k  j2,k 

problem  variables. 


B.U.F.  DESCRIPTION 


2d 


The  variables*  selected  to  describe  the  B.U.F.'s  per¬ 
formance  are : 

L  giving  the  number  of  B.U.F.**  involved  in  this  case 

[2  <  L  <  20]  . 

IBUD ***  indicating  the  nature  of  the  unloading  mode  at  the 

beach. 

IBCSL ****  indicating  which  of  the  B.U.F.  use  strategies  is  to 
be  used. 

TBC n  .  giving  the  number  of  units  of  time  after  the  start 
of  the  mission  that  the  Jlth  ( i  =•-  1,2... L)  B.U.F.  is 
expected  to  depart  from  its  base  [-999  <_  TBC ^  ^  <_ 
9999  (treated  as  a  floating  point  number)]. 


*The  magnitude  of  all  the  variables  selected  to  describe  the 
B.U.F.'s  performance  except  the  magnitude  of 

INBC  .  „  ( =  1,2.  .  .4,6,7.  ..  10;  l  =  1,2.  .  .1) 

is  affected  by  environment  state. 

**Each  B.U.F.  is  identified  by  a  distinct  number,  l,  such  that 
1  <_  l  <  L. 

***IBUD  -  1  if  the  unloading  mode  at  the  beach  is  to  be  in 
parallel. 

-2  if  the  unloading  mode  at  the  beach  is  to  be  se¬ 
quential. 

****If  IBCSL  -  1,  B.U.F.  use  strategy  SLBCA  is  used  to  select 

the  appropriate  B.U.F.  when  necessary. 

=  2,  B.U.F.  use  strategy  SLBCB  is  used  to  select 
the  appropriate  B.U.F.  when  necessary. 

=  3, 4...  9,  additional  B.U.F.  use  strategies  to  be 
developed  by  the  user  if  so  desired,  for 
selecting  the  appropriate  B.U.F.  when 
necessary. 
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TBC ^  ^  giving  the  time  required  for  the  Jlth  (l  =  1,2...L) 
3.U.F.  to  reach  point  A  on  the  beach  from  its  base 
[ 0 .  <  TBC-2  ^  £  9999  (treated,  as  a  floating  point 
number) ] . 

TBC _  .  giving  the  time  required  for  the  Jlth  ( l  =  1,2...L) 

•Jj  ) t 

B.U.F.  to  be  made  ready  to  start  the  unloading  opera¬ 
tion  after  the  Jlth  B.U.F.  has  arrived  at  point  A  on 
the  beach  [0.  <  TBC^  ^  <_  9999  (treated  as  a  floating 
*  point  number)]. 

TBC A  0  giving  the  time  required  for  the  Jlth  (l  =  1,2...  L ) 

B.U.F.  to  travel  to  the  Jlth  beach  unloading  area  from 
point  A  on  the  beach  after  the  Jlth  B.U.F.  has  been 
made  ready  to  travel  [0.  <  TBC 4  ^  <  9999  (treated  as 
a  floating  point  number)]. 

TBC g  ^  giving  the  time  required  to  release  any  of  the  pay- 
load  units  utilizing  the  Jlth  (i  =  1,3 ...L)  B.U.F. 

This  operation  is  performed  only  if  the  appropriate 
T.V.  is  beached  and  ready  to  commence  the  unloading 
operation  and  the  Jlth  B.U.F.  has  reached  the  Jlth 
beach  unloading  area.  If  that  is  not  the  case,  the 
releasing  of  the  payload  unit  must  wait  until  the  two 
requirements  given  above  are  satisfied  [ 0 .  <  TBCg  ^ 

<_  9999  (treated  as  a  floating  point  number)]. 

TBC „  »  giving  the  time  required  to  secure  any  of  the  payload 
unit  on  the  Jlth  (l  =  1,2...L)  B.U.F.  after  the  pay- 
load  unit  has  been  released  [0.  <  TBC?  ^  9999 

(treated  as  a  floating  point  number)]. 

TBC Q  .  giving  the  time  required  after  the  payload  unit  in 
question  has  been  secured  on  the  Jlth  (l  =  1,2...L) 
B.U.F.  to  (a)  transport  any  of  the  payload  units  from 
the  Jlth  beach  unloading  area  to  point  A  on  the  beach 
by  utilizing  the  Jlth  B.U.F.,  (b)  unload  and  free  the 


TBCa  o  payload  unit  from  the  ith  B.U.F.  and  (c)  make  the 
(cont'd.)  £th  B.U.F.  ready  to  travel  again  [ 0 .  <  TBCg  ^  <_ 

9939  (treated  as  a  floating  point  number)] . 

TBCa  .  giving  the  time  required  to.  prepare  the  ith 

(Si  =  i,2...L)  3.U.F.  for  departure  after  it  has 
terminated  its  mission  [ 0 .  <  TBC^  ^  <_  9999  (treated 
as  a  floating  point  number)]. 

TBCin  giving  the  time  required  for  the  £th  (Si  -  1,2...L) 

A.  U  j  X 

B.U.F.  to  reach  its  base  after  it  has  been  made 
ready  for  departure  [0 .  <  tbc2q  l  —  9999  (treated  as 
a  floating  point  number)] . 

IDBC .  *  indicating  the  nature  of  each  of  the  ft  (fit  = 

04  ,  *> 

1,2, 3, 4, 6,  7  , . .  .10)  processes  described  above  for 
each  of  the  l  (Si  =  1,2...  L)  B.U.F. 


INBC  .  0  giving  the  seed  of  the  U(031)  distribution  which 

serves  to  predict  the  stochastic  behavior  of  each  of 
the  ji*  (fn  =  1,2,3,  4,  6,7.  ..  10)  processes  described 
above  for  each  i  ( i  =  1,2...L),  if  IDBC ^  =  2 ** 

[1  <  INBC  .  .  <  10s  -1] . 

J1*  j  * 

AMINBC  .  .  giving  the  minimum  value  of  the  range  of  variation 

04j  “ 

of  TBC  .  .  (.jit  -  1,2, 3, 4, 6,7. .  .  1C;  Si  =  1,2...L),  if 


of  TBC 


04  t  & 


IDBC.  .  =  2**  [-999.99999  <  AMINBC ,  <  9999.99999 ; 

J  4  j  *  •*  J 

0.  <  AMINBC.  g  <  9999.99999,  j,  =  2, 3 , 4 , fi, 7 . . . 20 ] . 

03  j  * 

JfBC .  .  giving  the  maximum  value  of  the  range  of  variation 
04  J  * 

Of  TBC  ^  ^  (fit  =  1,2,  3,4 ,6 ,7 .  .  .10 ;  l  =  1,2... L ),  if 
IBBC^  2**  [-999.99999  <  AMAXBC  2  £  <  9999.99999; 

0.  <  AMINBC.  0  <  9999.99999,  js  =  2, 3,4,6,  7.  .  .  00]  . 

* IDBC  .  ,  »  2  if  the  process  is  deterministic. 

*  *  2  if  the  process  is  stochastic,  drawn  from 
0(0,1)  distribution. 

=  3,4. .. 9999  if  the  process  is  stochastic,  drawn 
from  distributions  to  be  developed  by  the  user, 
if  so  desired. 


**If  IDBC  jit  t  g=l,  then  INBCj 
are  not  problem  variables4/ 


AMINBC ,•  .  and  AMAXBC ,•  0 

i/lf  jX  i/4  J  X 


T.  V.  DESCRIPTION 


The  variables*  selected  to  describe  the  T.V. *s  perform¬ 
ance  ate: 

I  giving  the  number  of  T.V.**  involved  in  this  case 

(1  <  I  <  20). 

IWA1SL***  indicating  which  of  the  T.V.  use  strategies  is  to  be 
used  in  W.A.l. 

JVA25L****indicating  which  of  the  T.V.  use  strategies  is  to  be 
used  in  W.A.II. 

AMAXTV „  .  giving  the  weight  of  payload  together  with  that  of 

0  j  t  ^ 

the  associated  lashings  that  the  ith  (i  =  1,2... I) 
T.V.  can  carry  in  any  of  its  trips  [ 0 .  <  AMAXTV  . 

^  o  j  "V 

£ 9999  (treated  as  a  floating  point  number)}. 


*The  magnitude  of  the  variables  IWA1SL ,  IWA2SL,  AMAXTV ?  ^  (i  - 
1,2... I),  TTVj  i  (fa=l,2,  4.  .  .  7, 9, 10, 11, IS, 15, 16, 17),  • 

(37=1,2.  ..  5,9,10.  ..  13,15,16,1?)  AMINTVj,  and  AMAXTV j  ±  07  3 
is  affected  by  the  environment  state.  ’  The  magnituaVof  the 
variables  WT1MAX,  WT2MAX,  IDBRTVjSji  (39  =1,2. . .8)  BRKTVjg  i  is 
affected  by  the  environment  state  and  breakdown  considerations. 

**Each  T.V.  is  identified  by  a  distinct  number,  i,  such  that 
1  <  i  <  I. 

***If  IWA1SL=1,  T.V.  use  strategy  ASLTVA  is  used  to  select  the 

appropriate  T.V.  from  W.A.l  when  necessary. 

=2,  T.V.  use  strategy  ASLTVB  is  used  to  select  the 
appropriate  T.V.  from  W.A.l  when  necessary. 

=3, 4... 9,  additional  T.V.  use  strategies  to  be  de¬ 
veloped  by  the  user,  if  so  desired,  for  select¬ 
ing  the  appropriate  T.V.  from  W.A.l  when 
necessary. 

****If  IWA2SL=1,  T.V.  use  strategy  BSLTVA  is  used  to  select 

the  appropriate  T.V.  from  W.A.  II  when  neces¬ 
sary. 

« 2 ,  T.V,  use  strategy  BSLTVB  is  used  to  select 

the  appropriate  T.V.  from  W.A.  II  when  neces¬ 
sary. 

=3, 4...  9,  additional  T.V.  use  strategies  to  be 
developed  by  the  user,  if  so  desired,  for 
selecting  the  appropriate  T.V.  from  W.A.II 
when  necessary. 


AMAXTV „  .  giving  the  volume  of  payload  together  with  that  of 

•  J  ^ 

the  associated  lashings  that  the  ith  (i  =  1,2... I) 

T.V.  can  carry  in  any  of  its  trips.  Note  that  the 

definition  of  units  of  AMAXTV „  .  must  be  the  same  as 

k  ?>v 

that  of  VCn  (n  =  1,2. .kl2Nk)  [0  <  AMAXTV ?  i  <  99999 
(treated  as  a  floating  point  number)]. 

TTV  .  giving  the  number  of  units  of  time  after  the  start 

i  j  “t 

of  the  mission  that  the  ith  (i  =  1,2... I)  T.V.  is 
expected  to  depart  from  its  base  [-999  <_  TTV  ^  ^  < 
9999  (treated  as  a  floating  point  number) 


giving  the  time  required  for  the  ith  (i  =  1,2... I) 
T.V.  to  reach  W.A.I  from  its  base  [ 0 .  <  TTV n  .  < 
9999  (treated  as  a  floating  point  number)] 


giving  the  time  required  for  the  ith  (i  -  1,2... I) 
T.V.  to  reach  and  hook  up  to  any  of  the  ship  unload¬ 
ing  areas  from  W.A.I  and  prepare  the  ith  T.V.  for 
the  loading  operation.  Note  that  the  ith  T.V.  can 
leave  '^.A.I  only  when  there  is  a  ship  unloading  area 
free  to  receive  it.  [0.  <  TTV .  .  <  9999  (treated  as 

4,1, 

a  floating  point  number)]. 

giving  the  time  required  for  the  ith  (i  =  1,2... I) 
T.V.  to  refuel,  when  necessary,  after  it  has  hooked 

up  at  any  of  the  ship  unloading  areas  [0.  <  TTV ,  .  < 

0  }  1* 

9999  (treated  as  a  floating  point  number)]. 

giving  the  time  expected  for  the  ith  (i=l,  2. . .  I)  T.V. 
will  operate  at  other  than  zero  speed,  during  a 
complete  cycle.  [0.  <  TTV ^  ^  <_  9999  (treated  as  a 
floating  point  number)] . 


giving  the  time  expected  that  the  ith  (i  •  1,2... I) 
T.V.  will  operate  at,  other  than  zero  speed,  without 
refueling.  [0.  <  TTV „  .  <  9999  (treated  as  a  float- 
ing  number)] 


TTV  .  .  giving  the  time  required  for  any  payload  unit  to  be 

9  j  y>  ^ 

secured  on  the  ith  (i  =  1,2...!)  T.V.  after  it  has 
been  unloaded  into  the  ith  T.V.  and  freed  from  the 
appropriate  S.U.F.,  and  after  the  S.U.F.  in  question 
has  been  made  ready  to  travel  again  [G.  <  TTV.  .  < 

%7  j 

9999  (treated  as  a  floating  point  number)]. 

TTV,.  •  giving  the  time  required  for  the  ith  (i  =  1,2... I) 

T.V.  to  unhook  and  be  made  ready  to  travel  after  the 
last  payload  unit  unloaded  into  the  ith  T.V.  has  been 
properly  secured.  [0.  <_  TTV ^  ^  £  9999  (treated  as  a 
floating  point  number)]. 


TTV.  ,  .  giving  the  time  required  for  the  ith  (i  =  1,2... I) 

T.V.  to  reach  W.A.II  from  any  of  the  ship  unloading 
areas  [0 .  <  TTV..  .  <  9999  (treated  as  a  floating 

1  1  j  'l* 

point  number)]. 

TTV.,  .  giving  the  time  required  for  the  ith  (i  =  1,2... I) 

T.V.  to  reach  any  of  the  beach  unloading  areas  from 
W.A.II  and  then  beach  and  be  made  ready  for  the  un¬ 
loading  operation.  Note  that  the  ith  T.V.  can  leave 
W.A.II  only  when  there  is  a  beach  unloading  area  free 
to  receive  it  [0.  <  TTV.,  .  <  9999  (treated  as  a  float- 
ing  point  number)]. 

TTV .  c  .  giving  the  time  required  for  the  ith  (i  =  1,2... I) 

1  b  ,V 

T.V.  to  be  made  ready  to  travel  again  after  the  last 
payload  unit  carried  on  any  of  its  trips  has  been  se¬ 
cured  to  the  appropr  ate  B.U.F.  [0.  <  TTV..  „•  <  9099 

1  0  j  X-  — ■ 

(treated  as  a  floating  point  number)] . 

TTV..  .  giving  the  time  required  for  tht  ith  (i  =  1,2... I)  T.V. 

to  reach  W.A.I  from  any  of  the  beach  unloading  areas 

after  it  has  been  made  ready  to  travel  [0.  <  TTV..  .  < 

Iff,, i  - 

9999  (treated  as  a  floating  point  number)]. 


TTV 


17,  i 


IDTV 


INTV .  . 

J7  i1 


AMI NT  V 


J7 


AMAXTV. 


giving  the  time  required  for  the  ith  (i  *  1,2... I) 
T.V.  to  reach  its  base  from  any  of  the  beach  un¬ 
loading  areas  after  it  has  completed  its  mission 
[0 .  <  TTV .  „  .  <_  9999  (treated  as  a  floating  point 


number) ] . 

*  indicating  the  nature  of  each  of  the  j7  (1,2...  5,9, 
10 ...  12,15, 18,17 )  processes  described  above**  for 
each  of  the  i  (i  =  1,2...I)  T.V. 


giving  the  seed  of  the  U(0,1)  distribution  which 
serves  to  predict  the  stochastic  behavior  of  each 
of  the  J7  fj7  «  1,2. .  .5,9,10. .  .13,15,16,17)  pro¬ 
cesses  described  above  for  each  i  (i  -  1,2... I),  if 


IDTV 


=2***  [1 


<  INTV.  .  < 
«77  *  '*• 


10s  -1]  . 


giving  the  minimum  value  of  the  range  of  variation 


of  TTV  .  i  (j?  =  1,2..  .5,9,10.  .  .13,15,16,17;  i  - 
1,2...I)\  if  IDTV.  „.= 2***1-999.99999  <  AMIN  TV 


'7  J 


9999.99999;  0.  <  AMINTV .  . 

2.  .  .5,9,  10.  .  .  12,15,16,17]*. 


1,  i 

<  9999.99999,  j§  =  2, 


giving  the  maximum  value  of  the  range  of  variation 
(j7  =  1,2.  .  .5,9,10.  ..  12,15,16,17;  i  = 


of  TTV  .  . 

J7  j1 

1,2... I),  if  IDTV.  .  =  2***  [-999.99999  <  AMAXTV.  . 

67  !«• 

<  9999.99999,  j  = 


<  9999.  99999;  0.  <  AMAXTV ’ 


2,3.  .  .  5,  5,2(7.  .  .  12,  15,  16,17]  . 


*IDTV .  .-1  if  the  process  is  deterministic. 

=2  if  the  process  is  stochastic,  drawn  from  U(0,1) 
distribution. 

-2,4. . .9999  if  the  process  is  stochastic,  drawn  from 
distributions  to  be  developed  by  the  user,  if  so 
desired. 

**Process  3  is  the  process  associated  with  the  waiting  of  a 
T.V.  in  W.A.I. 

12  is  the  process  associated  with  the  waiting  of  a 
T.V.  in  W.A.II. 

***If  IDTV  .  >  ~  1 ,  then  INTV.  .,  AMINTV.  .  and  AMAXTV  .  . 

07>i  *  o7  >1  o7 

are  not  problem  variables. 


WT1MAX 


giving  the  maximum  time  that  any  of  the  j  T.V.  is 
expected  to  wait  in  W.A.I  at  any  time  during  the 
mission  [0.  <  WT1MAX  <  9999999. 99]  . 

WT2MAX  giving  the  maximum  time  that  any  of  the  j  T.V.  is 

expected  to  wait  in  W.A.II  at  any  time  during  the 
mission  [0 .  <  WT2MAX  <  9999999.99] . 

IDBRTV  .  .  *  indicating  the  presence  or  absence  of  breakdown 

J9  j1 

considerations  for  each  of  the  j3  (j9  =  1,2...  8) 
processes  (see  Table  3-1)  for  the  ith  (i  =  1,2. .  I) 
T.V.  and,  in  the  event  that  breakdown  considera¬ 
tions  are  present,  their  nature. 

INBRTV  .  .  giving  the  seed  of  the  U(0,1)  distribution  which 

J9  ,  ^ 

serves  to  predict  the  stochastic  behavior  of  each 
of  the  j9  (j9  -  1*2...  8)  processes  described 
above  for  each  i  (i  -  1,2...I),  if  IDBRTV .  2 ** 

[1  <  INBRTV .  .  <  109  -1].  9  J 

Jj  j1 

BRKTV .  .  giving  the  probability  that  a  breakdown  will  occur 

0  j  j1 

during  the  th  (jjo  =  1,2,4,6,7,81  process  for 
each  i  (i  =  1,2... I),  and  the  probability  that  a 
breakdown  will  occur  if  the  ith  T.V.  waited  WT1MAX 
or  more  units  of  time  in  W.A.I  during  the  2nd 
process***  and  the  probability  that  a  breakdown 
wilx  occur  if  the  ith  T.V.  waited  WT2MAX  or  more 
units  of  time  in  W.A.II  during  the  5th  process***, 
if  IDBRTV.  .  =2**  [ 0 .  <  BRKTV.  .  <  1.0000 ]. 

1/9  ,  i  J 9  >  i 

*  IDBRTV .  .  «»  1  if  there  are  no  breakdown  considerations . 

^9  *  2  if  there  are  breakdown  considerations  which 

are  drawn  from  a  U(0,1 )  distribution. 

=  2, 4...  99  if  there  are  breakdown  considerations 

which  are  drawn  from  distributions  to  be  devel¬ 
oped  by  the  user,  if  so  desired. 

**  If  IDBRTV ,i~l,  t'lienlNBRTV j  ^  and  BRKTV j9  si  are  not  prob¬ 
lem  variables. 

***  If  the  T.V.  waited  less  the  probability  is  scaled  down 
linearly. 


the  breakdown  considerations  of  each  T.V.  regarding 


oo  eu 

0  o. 

43  -h  oo 
•H  U  -H 
Vl  VI  A 
U  VI 

at  oo 

V  U  IH 
Q  44  O 


0) 

. 

a 

• 

3 

• 

c  • 

a 

C  TO 

3  TO 

o 

3  to 

® 

03 

3  V 

co  o 

a) 

03 

0  TO 

0) 

CD 

O  0 

C‘ 

o  > 

XJ  T3 

> 

XJ 

CO 

VI  M 

> 

XJ 

TO  > 

o 

•H 

TO  r4 

4-1  H 

•H 

CO 

4  OH 

i4  40  H 

Vi 

^  O 

10 

o 

a) 

0 

o 

0  O 

eo  a. 

00 

<0  > 

60  V 

> 

00 

a 

60  Vl 

> 

60 

0  > 

c 

d 

®  c 

a  Vi 

c 

c 

o 

C  43 

C 

C 

n  G 

•H  oo 

•H 

Vi  -H 

•H  A 

•H 

•H 

u 

44 

•H 

•r4 

A  44 

■a  -h 

43 

TO 

a. 

TO  0 

TO 

c  x: 

• 

H 

• 

Vl  0 

• 

M 

• 

Vl  A 

• 

Vl 

0  • 

co  vi 

0) 

c0 

0)  > 

0  .G  > 

c0 

CD  0) 

0  H  > 

0  43  > 

60 

0 

60  £  • 

00  H 

• 

60  t4  0 

60 

• 

60  H  • 

0  4-1 

•H 

a) 

H  H 

0 

H 

a) 

40  -H 

0 

H 

a> 

R 

M  O 

VI 

u 

Vl 

VI  VI 

Vl  • 

u 

<U 

• 

O 

1-4 

a) 

.  4) 

•  00 

60 

• 

•  JC 

•  a 

J1 

• 

44  60 

•  W 

ji 

• 

0  A 

>  C 

C 

> 

M  XJ 

>  <0 

xj 

> 

o  a 

>  • 

XJ 

> 

0  Vi 

•  o 

•H 

• 

• 

•  0 

• 

44 

•  ■< 

• 

0 

H  -H 

VI 

H  <0  4-1 

H  U 

CH 

H 

0  VI 

H  • 

CH 

H 

Vl  44 

VI 

•H 

•  o 

4 

O 

fi  44 

3: 

O 

0  J 

£  <0 

00 

A  3 

X 

4= 

O  4 

A 

x: 

U  Vc 

3 

a 

c 

CJ  60 

d 

a 

44  3 

U  0 

d 

a 

60  G 

«  <U 

cd 

0  O 

4  C 

o 

co 

VI 

0  O 

o 

cd 

G  O 

<0  -a 

0 

a) 

O  -H 

0  -H 

•H 

a) 

0  0 

0  Vl 

a) 

•H  -t4 

•H 

A 

VI  VI 

TO 

U 

VI  A 

44 

XJ 

TO  VI 

4-1  00 

VI 

M-l 

4-1  O 

44  0 

a 

m 

0  VI 

44 

a 

M-l 

0  U 

O  C 

O 

c 

o  o 

d 

o 

TO 

O  0 

d 

o 

O  G 

O 

•a 

oo  3 

r4 

3 

•H  TO 

0 

d 

i-4  3 

00  O 

c 

C/3 

10  4-1 

0  C 

4-1 

CO 

0  a 

0  0 

44 

CD 

G  41 

c 

■0 

d 

a> 

C  3 

c 

C  4 

a  u 

d 

3 

o  c 

o 

Vi  co 

O 

co 

o 

o 

O  0 

c0 

o 

0 

•h  s  to 

•H 

as 

ti  a 

•H 

O  TO 

•H 

•H 

A 

vi  o 

0 

XJ 

a) 

VI  -r4 

a) 

XJ 

0 

vl  60 

OJ 

XJ 

U  0 

<0  TO 

> 

d 

60  Vi 

4  jC 

u 

co 

c  > 

0  e 

u 

CO 

0  Vl 

U  44  .-4 

U 

c  to 

W  0 

C0 

u 

3  iH 

Vl  *H 

cO 

U 

0  uS 

0  CO 

o 

a) 

•H 

0 

0) 

o  o 

0  TO 

a) 

43 

TO  0 

> 

•a 

TO  0) 

TO  0 

03 

T3  > 

TO  0 

03 

0 

•H  Vi 

c 

•H 

io  oo 

•H  A 

C/3 

•H 

a 

•H  O 

CD 

•H 

0  0 

00  A 

•H 

03 

O  0) 

0  VI 

0) 

CD 

0  44 

0  1-4 

a) 

03 

A  0 

C 

C 

.H  CJ 

a 

c 

0 

a  c 

u 

d 

VI  u 

O  0) 

• 

o 

C  O 

O  44 

o 

o 

u  • 

O  3 

o 

o 

o 

o  A 

> 

o 

3  C 

U  O 

V4 

u 

43  > 

o 

u 

a 

44  Vl 

H 

• 

a 

a. 

• 

A 

a. 

O  0. 

C 

H 

c 

a 

a  ►. 

c 

0)  H 

a  u 

c 

3 

3  n  co 

3  C 

CD 

3  A 

3  0 

03 

^4  0 

o  • 

0 

o 

A  -H 

O  4 

•H 

o 

H  0 

O  0 

•H 

o 

G  44 

•a  m 

£ 

•a 

00  A 

TO 

X 

•a 

A 

TO  43  A 

0  A 

x  • 

VI 

A 

VI 

4*  0 

4J 

VI 

VI 

VI 

00  4jJ 

to 

V 

4  O 

cO 

• 

0  0 

c0 

0 

0  • 

44 

0 

A  44 

0  Vl 

CW 

a) 

W  44 

0  A 

44 

a> 

O  41 

VI  3 

O 

Vi 

VI  O 

Vl  44 

O 

Wi 

H  O 

Vl  VI 

o 

u 

u  o 

43 

A 

A 

X) 

• 

A 

XJ 

41 

c 

c 

44  61 

H 

03 

<  c 

41 

0 

0 

0  -H 

o 

0 

o  c 

0  H 

C 

0/ 

•  o 

0  O 

G 

a> 

H  G 

A 

•H 

A 

o 

A  • 

o 

£344 

A 

o 

•  O 

xj  oo  xj 

VI 

►.  v4 

VI  < 

•H 

XJ 

VI 

vi 

i4 

V)  -<  -H 

a 

a 

C  VI 

• 

Xi 

c  u 

G 

VI 

•  XJ 

01  *rl 

a 

<0 

co  CO 

0  3 

co 

(0 

44  C 

0  0 

0 

0  3  0 

0  VI 

3 

V 

Vi 

0 

u 

0) 

3 

0 

Vl 

a> 

u 

A  *H 

M-l 

A 

O  V 

43  O 

a) 

XJ 

60  41 

XI  o 

0 

XJ 

O  0) 

•H  10 

•H 

VI  TO 

•H  VI 

•a 

*H 

C 

44  VI 

TO 

•H 

XJ  T3 

vi  3 

c0 

Vi 

t4 

Vl 

•H 

44  0 

Vl 

i4 

•H 

o 

U 

0.  61 

O  0. 

CO 

a 

VI 

u  a. 

0 

a 

a.  cd 

0)  00 

0) 

61 

•H  C 

0  44 

d 

co 

•H  0 

0  44 

G 

CD 

•h  d 

a>  vi 

M 

V 

Vi  O 

0  Vl 

0 

a) 

0  Vl 

0  Vl 

o 

a> 

u  o 

Q  4-1 

cO 

a 

vi  CJ 

Q  vi 

a 

Q 

3  0 

a  « 

CJ 

Q 

xj  a 

CM 

CO 

Co 

Cs. 

II 

II 

n 

II 

II 

II 

c 

3 

O 

•o  • 
M  -a 
0  0  0 
VI  «  > 
•H  M  r4 
A  O 
00  > 

a  ®  c 

•H  JS  *H 
TO  H 
M 

0  > 

60  •  • 

0  oo  r 

M  0 

V  V 
•  H  A 

>  *  « 

H  00  4-i 

C  O 
A  -H 
CJ  TO  G 
0  0  O 
0)  O  •H 
r-C  VI 
n  C  U 

o  n  a 

3 

a)  40  n 

g  o 

O  «  10 
-H  0 
4J  43  0) 
CO  Vi 
Vi  O  O 

a>  a 
to  *i  co 

•H  00 

mm  v 
coo 
o  o 
o  >,  c 
a  a 

C  <0 
3  00 
O  0  -H 
O  0£ 
^  V  U 
co  44 
a)  44 
coo 
A  6) 

CO  (0 

0  xi  C 

A  o 

4J  00  "H 

VI  VI 
00  -H  CO 
0  U 
A  O  0 
•H  Ufl 
Vi  -H 
U  0.00 

CO  <H  c 

«  M  O 

a  vi  o 


With  the  above,  the  list  of  the  variables  that  are 
utilized  to  define  the  five  subsystems  involved  in  our  study 
and  at  the  same  time  appear  as  input  in  our  computer  program 
is  complete.  There  exist  two  additional  inputs  to  the  com¬ 
puter  program  which  serve  to  control  the  program's  performance 
and  which  for  the  sake  of  completeness  we  include  here.  These 
are: 

NCASES  giving  the  number  of  cases  to  be  processed  in 
each  program  execution  [1  <_  NCASES  _<  99],  and 
NRUNS  giving  the  number  of  runs  to  be  processed  for 

the  j^th  (j  =  1,2. . .NCASES)  case  [1  <  NRUN  £  100] 

We  continue  now  by  listing  the  remaining  variables  neces¬ 
sary  to  complete  each  subsystem's  description. 

MOTHER  SHIP  DESCRIPTION 

The  additional  variables  selected  to  complete  the  descrip 
tion  of  the  Mother  Ship's  performance  are: 

TAMP ,  TIP ,  T2P*  giving  the  time  fluctuation  associated  with 
TAM,  T1  and  T2  respectively. 

TTM  giving  the  total  time  taken  by  the  Mother  Ship  to  com¬ 
plete  its  operation.  This  includes  the  time  taken  to 
free  the  Mother  Ship  from  its  mooring  and  to  make  it 
ready  for  travel  again  after  all  S.U.F.  are  secured  in 
position  and  all  T.V.  have  cleared  the  Mother  Ship. 

If  IDMA  and/or  ID1  and/or  ID2  equal  1 ,  then  TAMP  and/or  TIP 
and/or  T2P  equal  zero  respectively.  If  that  is  not  the  case, 
TAMP ,  TIP  and  T2P  are  determined  by  drawing  from  the  appro¬ 
priate  distribution,  as  dictated  by  the  magnitude  of  IDMA, 

ID1  and  ID2  respectively. 


PA'ILOAD  DESCRIPTION 


The  variables  given  above  suffice  to  describe  the  pay- 
load  characteristics  to  the  degree  of  accuracy  and  extent 
desired ,  and  so  no  additional  variables  are  needed. 

S.U.F,  DESCRIPTION 

The  additional  variables  selected  to  complete  the  descrip 
tion  of  the  S.U.F.'s  performance  are: 

TSC6  .  giving  the  time  that  the  fcth  (k  =  1,2... K)  S.U.F. 

has  to  wait  in  the  kth  ship  unloading  area  before 
it  can  unload  the  payload  unit  that  it  is  trans¬ 
porting.  The  kth  S.U.F.  has  to  wait  until  the 
appropriate  T.V.  is  properly  secured  in  the  fcth 
ship  unloading  area  or  until  the  previously  un¬ 
loaded  payload  unit  is  properly  secured  in  the  T.V. 
in  question. 

TSCP.  .*  giving  the  time  fluctuation  associated  with  TSC  .  , 

*7  2  •*  *  <7,  » 

(j  2  =  1,2.  .  .5,7,  8,  9;  k  -  1,2...K).  ■ 

AMAXSCq  k  giving  the  total  time  taken  by  the  fcth  (k  -  1,2. . ,K) 
S.U.F.  to  complete  its  mission.  This  includes  the 
time  taken  to  secure  the  S.U.F.  in  its  original 
position. 


*If  IDSC.  .  equals  1,  then  TSCP.  ,  equals  zero.  If  that  is 
<7  2  *  K  3z 

not  the  case,  TSCP.  ,  is  determined  by  drawing  from  the 

Q2  ,  K 

appropriate  distribution,  as  dictated  by  the  magnitude  of 


IDSC  . 

Ci  » 


km 


B.U.F.  DESCRIPTION 


The  additional  variables  selected  to  complete  the 
description  of  the  B.U.F.'s  performance  are: 

TBC ^  ^  giving  the  time  that  the  Jlth  (l  =  i,2...L )  B.U.F. 

has  to  wait  in  any  of  the  beach  unloading  areas 
before  it  can  release  the  appropriate  payload 
unit.  The  Jlth  B.U.F.  has  to  wait  in  a  beach  un¬ 
loading  area  until  the  appropriate  T.V.  is  properly 
beached  and  has  been  made  ready  for  the  unloading 
opeiation. 


TBCP .  0*  giving  the  time  fluctuation  associated  with  TBC  .  n 

i  x  j  * 

(j\  «  1,2,3,4,6,7.  .  .  10;  l  -  1,2. ..L). 


AMAXBC 5  ^  giving  the  total  time  taken  by  the  Jlth  (i.  =  1,2...L) 
B.U.F.  to  complete  its  mission.  This  includes  the 
time  taken  for  the  Jlth  B.U.F.  to  reach  its  base. 


*If  IDBC  .  „  equals  1,  then  TBCP.  0  equals  zero.  If  that  is 

not  the  case,  TBCP.  «  is  determined  by  drawing  from  the 

<7  4  i  * 

appropriate  distribution,  as  dictated  by  the  magnitude  of 


IDBC  . 


Jl* 


T .  V.  DESCRIPTION 


is  4 


The  additional  variables  selected  to  complete  the  de¬ 
scription  of  the  T.V. 1 s  performance  are: 

TTV 7  .  giving  the  time  that  the  ith  (i  -  1,2... I)  T.V. 

has  to  wait  in  W.A.I.  The  ith  T.V.  has  to  wait  in 
W.A.I  until  the  appropriate*  ship  unloading  area 
becomes  free  to  receive  it.  At  the  start  of  the 
mission/  as  soon  as  the  Mother  Ship  is  moored/  all 
ship  unloading  areas  become  free.  Subsequently/  a 
ship  unloading  area  becomes  free  as  soon  as  the 
T.V.  that  is  being  served  alongside  the  ship  unload¬ 
ing  area  in  question  is  unhooked  and  made  ready  to 
commence  its  journey  to  W.A.II. 

TTVg  £  giving  the  time  that  the  ith  (i  -  1,2... I)  T.V.  has 
to  wait  in  any  of  the  ship  unloading  areas  awaiting 
the  appropriate  S.U.F.'s  arrival. 

TTV  0  .  giving  the  time  that  the  ith  (i  -  1,2... I)  T.V.  has 

1  u  j 

to  wait  in  W.A.II.  The  ith  T.V.  has  to  wait  in 
W.A.II  until  the  appropriate**  beach  unloading  area 
becomes  free  to  receive  it.  At  the  start  of  the 
mission/  as  soon  as  the  Jlth  ( SL  =  1,2...L)  B.U.F. 
arrives  at  point  A ,  the  ith  beach  unloading  area  be¬ 
comes  free.  Subsequently  a  beach  unloading  area 
becomes  free  as  soon  as  the  T.V.  that  is  being 
served  at  the  beach  unloading  area  in  question  is 
made  ready  to  travel  again  for  W.A.I  or  its  base. 

*  Note  that  during  the  sequential  loading  mode  if  one  ship  un¬ 
loading  area  is  not  free,  then  all  ship  unloading  areas  are 
considered  busy. 

**Note  that  during  the  sequential  unloading  mode  if  one  beach 
unloading  area  is  not  free  then  all  beach  unloading  areas 
are  considered  busy. 


TTV 14  ^  giving  the  time  that  the  ith  (i  =  1,2... I)  T.V. 

has  to  wait  in  any  of  the  beach  unloading  areas 
awaiting  the  appropriate  B.U.F.'s  arrival. 

TTVP .  . *  giving  the  time  fluctuation  associated  with  TTV  . 

(j7  =  1,2.  .  .5,9,10.  .  .'13,15,16,17;  i  =  1,2.  ..I).  7  * 

AMAXTV a  .  giving  the  total  time  taken  by  the  ith  (i  = 1,2. . .1) 

O  j  "V 

T.V.  to  complete  its  mission.  This  includes  the 
time  taken  for  the  ith  T.V.  to  reach  its  base. 

The  above  completes  the  list  of  variables  that  we  se¬ 
lected  to  describe  our  system  for  this  study.  In  the  contin¬ 
uation  of  this  study,  the  list  will  be  augmented  as  necessary 
to  deal  with  the  statistical  considerations  and  the  concept  of 
antithetic  variance.  These  considerations  and  this  concept 
will  be  presented  in  a  separate  report  which  will  follow  the 
present  one.  The  above  completes  the  detailed  discussion  of 
the  first  topic  of  the  Problem  Definition.  We  now  proceed 
with  the  second  topic  of  the  Problem  Definition,  which  deals 
with  the  establishment  of  the  relations  among  the  problem 
variables . 


*If  IDTV .  .  equals  1,  then  TTVP.  .  equals  zero.  If  that  is 

3 7  3y  »  ^ 

not  the  case,  TTVP .  .  is  determined  by  drawing  from  the 

07  »  ^ 

appropriate  distribution,  as  dictated  by  the  magnitude  of 
IDTV .  .. 


ii)  Establishment  of  the  Relations  Among  the 
Problem  Variables 

Because  in  the  fourth  section  we  will  undertake  to  con 
struct  the  mathematical  model  of  our  system,  we  should,  be¬ 
fore  we  proceed,  establish  all  the  relationships  that  exist 
among  the  variables  that  we  selected  to  describe  our  system 
This  is  because  the  mathematical  model  is,  by  definition, 
the  replica  of  our  system  in  the  form  of  mathematical  equa¬ 
tions,  and  therefore  any  such  relations  must  be  part  of  it 
in  order  to  permit  it  to  be  a  true  replica  of  the  original. 
So  we  proceed  by  establishing  all  such  relationships. 

(3.: 

TAMP  =  0  if  IDMA  =  1 

*  ( AMAXM-AMIRM) AR  +  AMINM  if  IDMA  =  2 

where  R  is  a  random  number  generated  from  a  U(031) 
distribution  utilizing  I RMS  as  the  first  seed  and 
then  the  updated  INMS  in  subsequent  generations. 

TIP  =  0  if  ID1  =  1  (3.; 

*  (AMAX1- AMIRI )*R  +  AMIRI  if  ID1  *  2 

where  R  is  a  random  number  generated  from  a  V(Q,1) 
distribution  utilizing  IR1  as  the  first  seed  and 
then  the  updated  IR1  in  subsequent  generations. 

T2P  =  0  if  ID2  =  1  (3.: 

=  (AMAX2-AMIR2) *R  +  AMIR2  if  ID2  =  2 

where  R  is  a  random  number  generated  from  a  U(031) 
distribution  utilizing  IR2  as  the  first  seed  and 
then  the  updated  IR2  in  subsequent  generations. 


TTM  =  TAM  +  TAMP  +  T1  +  TIP  +  T2  +  T2P  + 


max  < 


max  f AMAXSC 

if  (  £ 


max  [TSCP6'k] 


*  * 


(3.4) 


where  TSCP6  ^  records  the  last  time  that  the 
kth  (k  =  1,2... K)  ship  unloading  area  became  free. 


TSCP .  m  o  if  IDSC  .  !=  1 

J2*k  J2>k 


(3.5) 


=  [AMAXSC .  ,  -  AMINS C 


l 


j 2,  k 


'*•*)  4  * + 


AMINSC  .  - 

J2,k 


if  IDSC.  .  =  2 

C  23< 

where  J?  is  a  random  number  generated  from  a  U(0,1) 

distribution  utilizing  I'SCj  ka2  =  1,2...  5,  7,  8,  9; 

k  -  1,2. . . K)  as  the  first  seed  and  then  the  updated 

INSC.  .  in  subsequent  generations. 

j  2,k 


TSC6  £ ,  AMAXSC 6  k  The  relationships  of  TSC6  AMAXSC 6  k  (3.6, 

with  the  other  problem  variables  are  too  cumber- 
some  to  be  written  in  the  form  of  mathematical  equa¬ 
tions,  and  so  the  reader  is  referred  to  the  listing  of 
the  computer  program,  where  the  relationships  in  ques¬ 
tion  are  given  in  the  form  of  computer  coding.  The 
mathematical  form  of  these  relationships  can  be  ob¬ 
tained  from  the  computer  coding,  if  so  desired,  and  it 
is  similar  to  the  form  of  equation  (3.4). 

* 

The  magnitude  of  AMAXSCg  y  at  the  conclusion  of  the  mission 
of  the  kth  S.U.F.  J 

**The  magnitude  of  TSCPg  %  immediately  after  the  last  of  the 
N^  payload  units  has  Jbeen  unloaded  and  secured  in  the 
appropriate  T.V.,  and  after  the  T.V.  in  question  has  been 
unhooked  and  made  ready  to  depart  for  W.A.II. 


TBCPj  o 
0  *• 


0  if  IDBC  .  .  =  1 

o4A 


(3.8) 


=  \AMAXBC .  .  -  AMINBC .  .)  *  R  +  AMINBC  .  . 

[  '}4i*'  J4jX,J  J4'X’ 

if  IDBC;  o  =  2 

where  R  is  a  random  number  generated  from  a 
0(0,1)  distribution  utilizing  INBC  •  0 

U  4»  *• 

(04  -  1,2, 3,4, 6,7 ...  10;  l  «  1,2...L)  as  the 

first  seed  and  then  the  updated  INBC .  , 

0  4* 

in  subsequent  generations. 


TBC5  v  AMAXBCS  a  (3.9,10) 

The  comments  given  for  equations  (3.6)  and 
(3.7)  apply  also  for  equations  (3.9)  and  (3.10). 


TTVP  .  .  «  0  if  IDTV .  .  =  1  (3.11) 

0  1*  0  y  I  T* 

=  f AMAXTV •  •  -  AMINTV  •  ■ )  *  R  +  AMINTV •  ■ 

if  IDTV  ■  •  -  2 

d  J  J  ^ 

where  I  is  a  random  number  generated  from  a 
0(0,1)  distribution  utilizing  INTV^^i 
(3?  -  1,2.  .  .5,9,10.  .  .13,15,16,17;  i  =  1,2.  ..I) 
as  the  first  seed  and  then  the  updated  INTV 
in  subsequent  generations. 

TTVj  .,  AMAXTV 0  .  (3.12-16) 

0  2 1  •* ^  0 ,% 

j21  =  3,8,12,14 

The  comments  given  for  equations  (3.6)  and 
(3.7)  apply  also  for  equations  (3.12) -(3.16) 


% 


The  above  equations  establish  all.  the  relationships 
that  exist  among  the  problem  variables.  Now  we  proceed  with 
the  discussion  of  the  third  topic  in  the  Problem  Definition, 
namely,  the  identification  of  the  dependent  and  independent 
variables . 


iii)  Identification  of  the  Dependent  and 

Independent  Variables 

For  reasons  given  in  the  previous  section  it  is 
necessary,  before  we  proceed  any  further, -to  classify  the 
problem  variables  into  dependent  and  independent  ones.  As  a 
reminder,  the  independent  variables  are  those  variables  which 
must  be  prescribed  to  completely  describe  our  system  as 
desired.  The  dependent  variables  are  the  ones  remaining  in 
our  original  list  of  variables  after  the  independent  ones  have 
been  removed.  The  value  of  the  dependent  variables  can  be  ob¬ 
tained  from  the  variable  interrelationships. 

The  variable  classification  of  dependent  and  in¬ 
dependent  variables  can  be  accomplished  in  the  following  man¬ 
ner: 

For  each  interrelation,  one  of  the  variables  involved  is 
classed  as  a  dependent  variable  and  all  the  others  are 
classed  as  independent  variables.  With  this  method,  if 
the  same  variable  is  involved  in  more  than  one  inter¬ 
relation,  it  cannot  be  selected  to  serve  as  a  dependent 
variable  more  than  once. 


In  theory,  any  of  the  variables  involved  in  an  inter¬ 
relation  can  be  classed  as  a  dependent  variable.  However,  in 
practice,  whenever  it  is  possible  we  usually  attempt  to  class 
as  dependent  variables  the  ones  that  allow  us  to  simplify  our 
methodology  as  much  as  possible,  and  we  attempt  to  class  as 


I 


independent  variables  the  ones  which  are  the  most  meaningful 
to  the  user  so  he  can  assign  their  values  with  confidence. 
Whenever,  as  in  our  case,  that  is  impossible  because  of  the 
nature  of  the  interrelations,  we  usually  must  compromise. 

The  nature  of  the  equations  (3.1) -(3. 16)  forces  us  to  class 
the  variables  appearing  on  the  left-hand  side  of  these  equa¬ 
tions  as  the  dependent  variations.  However,  although  we  have 
no  choice  in  our  selection,  there  was  no  compromise  at  all,  as 
the  resulting  independent  variables  of  our  problem  are  both 
the  most  meaningful  of  the  variables  to  the  user,  and  at  the 
same  time  they  allow  us  to  simplify  our  methodology  the  most. 


So  our  dependent  variables  are: 


TAMP,  TIP ,  T2P,  TTM , 


TSCP  . 


TBCP  . 
*4* 

TTVP  . 

3 


k • 
J V 
V 


TSC6  k ,  AMAXSC6  k  (3 2=1,2.  .  .  5,7 ,8,9;k=l,2.  .  .K) 

TBC 5  ^  AMAXBC5  ^  (j 4=1,  2,3,  4,6,  7 .  .  .  10  ;i=l,2.  .  .  L) 
TTV  .  .,  AMAXTV Q  ,  (j „=1, 2.  .  .  5,9.10.  ,  .  13,15,16,17; 

3  j  j ,  *  0,1  / 

0 11  =  3,3,12,14;  i=l,2.  .  .1) 


and  our  independent  variables  are: 


TAM,  Tl,  T2,  IDMA,  ID1,  ID2,  ISMS,  AKINM,  AMAXM ,  INI,  AMIN1, 
AMAX1,  IN 2,  AMIN 2,  AMAX2, 

Nk,  WCn,  VCn  (k-1,  2.  .  .  K;  n=l,2...kl1  Nk) 


K, 

ISUD,  ISCSL, 

ISC3  2>k*  IDSC3  2>ki 

INSCo2,k >  AMINSCj2,k> 

AMAXSCj  i , 

3  2*k 

( 3  2=  1,2.  ..5,7,  3,9; 

k=l,2. . .K) 

2, 

IBUD ,  IBCSL, 

TBC3  4,1*  IDB C 3  4 ,  £•» 

INBCU,%,  AMINBCj 4}  i. 

AMAXBCj4}  j i 

( 3  £  =  1,2,  3,  4,  6,  7.. 

.10;  1=1,2. ..L) 

I.  IWA1SL ,  IWA2SL,  AMAXTV e  .,  AMAXTV ,  TTV  .  .,  IDTV  .  . 

6ji  7,i  3q>'1  d-psi 

INTV.  AMIN  TV .  AMAXTV.  .,  WT1MAX,  WT2MAX, 

3  7 «7  1  3  •pi'1’ 

IDBHTV .  .,  INBRTV  .  .,  BRKTV  .  . 

3  g*l  3  gj  ^  3  g  ,  I* 

(3  6=1,  2,  4,  5,6,  7,  9,  10,  11,12,15,16,17;  j?=l,2.  ..5,9,10... 
. . . 12,15,16,17 i  j g  —  1 , 2 ...  8 i  i=l, 2 ...  I) 

TIME,  WGHT  and  VOL. 


The  above  completes  our  discussion  about  the  identi¬ 
fication  of  the  dependent  and  independent  variables  of  our 
problem,  and  now  we  may  proceed  with  discussion  of  the  estab¬ 
lishment  of  the  range  of  variation  of  our  problem's  parameters. 


iv)  Establishment  of  the  Range  of  Variation  of  the 
Problem  Parameters 

For  the  reasons  given  in  the  previous  section,  it  was 
found  necessary  to  introduce  limitations  on  the  range  of 
variation  of  our  problem's  parameters.  The  limitations  intro¬ 
duced  in  our  study*  are  enumerated  and  explained  below. 

Methods  for  alleviating  each  limitation  when  it  is  found  un¬ 
acceptable  are  also  given  below,  whenever  it  is  deemed  neces¬ 
sary.  For  clearer  understanding  of  certain  of  the  limitations 
listed  below,  the  user  is  advised  to  refer  to  Fig.  3.1,  where 
the  subsys terns  involved  in  our  study  are  shown  diagrammati- 
cally . 

1.  The  entire  operation  is  assumed  to  have  started**  either 
when  the  mother  ship  arrives  at  the  theater  of  operations 
or  when  one  or  more  T.V.  or  B.U.F.  start  from  their  bases, 
whichever  is  sooner.  This  limitation  can  be  alleviated 
easily  by  assuming  that  the  entire  operation  commences  at 
any  prespecified  time,  as  desired. 

2.  The  entire  operation  finishes  when  the  mother  ship  is  freed 
from  its  moorings  and  is  ready  to  depart,  and  all  the  T.V. 
and  B.U.F.  reach  their  appropriate  bases.  This  limitation 


*Please  note  that  great  care  was  ^nken  in  introducing  the 
limitations  on  the  range  of  variation  of  our  problem's 
parameters  so  that  they  satisfy  all  the  requirements  about 
model  universality,  etc.,  given  in  Section  2  of  the  report. 

**At  the  start  of  the  operation  the  time  counter  in  our  pro¬ 
gram  is  initialized  to  zero. 


can  be  alleviated  easily  by  assuming  that  the  entire 
operation  terminates  at  any  prespecified  time,  as  desired. 

3.  As  was  already  noted  in  the  discussion  of  the  first  topic 

of  this  section,  limitations  of  the  form  u  .  <  u  < u 

amin  —  v  —  smax 

(where  y  is  any  of  our  independent  variables  or  NCASES  and 
NRUNS)  were  introduced  controlling  the  magnitude  variation 
of  the  different  problem  parameters.  These  limitations  can 
be  alleviated  by  changing  the  format  and/or  the  dimension 
statements  in  the  computer  program. 

4.  The  mission  of  each  S.U.F.  is  such  that  it  does  not  inter¬ 
fere  with  that  of  any  other  S.U.F. 

5.  Each  of  the  K  S.U.F.  is  assumed  to  complete  its  mission 
without  any  technical  difficulties.  This  limitation  can  be 
alleviated  by  introducing  breakdown  considerations  similar 
to  those  introduced  for  the  T.V. 

6.  Each  of  the  K  S.U.F.  requires  no  refuelling  during  its 
entire  operation. 

7.  The  unloading  at  the  ship  may  be  sequential  or  in  parallel, 
but  not  both.  In  addition,  the  unloading  mode  at  the  ship 
is  independent  of  that  at  the  beach.  Moreover,  in  the 
parallel  unloading  mode,  the  S.U.F.  performance  is  not  in¬ 
fluenced  by  the  number  of  S.U.F.  that  happen  to  be  unload¬ 
ing  at  any  instant  of  time. 


8.  The  tiroes  TSC0  ,  .  ,  .  depend  on  the  S.U.F.  characteris- 
tics  only,  and  are  independent  of  the  payload  unit  tney 
are  to  service. 

9.  The  time  TSC 7  ^  depends  on  the  S.U.F.  characteristics  only, 
and  is  independent  of  the  payload  unit  that  is  being  un¬ 
loaded  and  the  T.V.  into  which  the  payload  unit  in  ques¬ 
tion  is  being  loaded. 

10.  The  payload  units  capable  of  self -induced  rolling  motion 

may  be  loaded  into  the  T.V.  by  ramps.  When  this  is  done, 
the  ramps  are  treated  as  regular  S.U.F.  and  the  times 
TSC  .  =  2,  2.  .  .  £  ,7 ,  8 ,9 ;  k'-I.D.  of  ramps)  of  each 
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ramp  are  determined  by  averaging  the  unloading  character¬ 
istics  of  the  payload  units  that  are  to  be  unloaded  by  the 
ramp  in  question. 

11.  Each  payload  unit  loaded  at  the  area  associated  with  the 
kth  (k  =  1,2... K)  S.U.F.  can  be  unloaded  only  by  the  kth 
S.U.F. 

12.  The  kth  (k  =  1,2... K)  S.U.F.  is  assumed  to  be  able  to 
handle  the  heaviest  and  bulkiest  payload  unit  that  is  load¬ 
ed  in  the  area  that  is  associated  with  the  S.U.F.  in 
question. 

13.  Only  one  payload  unit  may  be  handled  by  a  S.U.F.  at  any 


one  time. 


14.  Once  a  T.V.  nooks  up  at  the  kth  (k  =  lt2...K)  ship  unload 
ing  area#  only  the  kth  S.U.F.  may  load  it. 

15.  Each  of  the  I  T.V.  may  hook  up  at  one  ship  unloading  area 
only  during  the  execution  of  any  one  of  its  trips . 

16.  Each  of  the  K  S.U.F.  continues  to  load  in  any  of  the  I 
T.V.  during  the  execution  of  any  one  of  its  trips  until 

i)  no  more  payload  units  can  be  loaded  into  the  T.V. 

in  question/  because  the  weight  capacity  of  the  T.V. 
will  be  exceeded,  or 

ii)  no  more  payload  units  can  be  loaded  into  the  T.V. 
in  question  because  the  volume  capacity  of  the  T.V. 
will  be  exceeded,  or 

iii)  there  are  no  more  payload  units  to  be  unloaded  by  the 
S.U.F.  in  question,  or 

iv)  any  combination  of  the  above  statements  becomes  true. 

17.  The  mission  of  t'zch  B.U.F.  is  such  that  it  does  not  inter 
fere  with  that  of  any  other  B.U.F. 

18.  Each  of  the  L  B.U.F.  is  assumed  to  complete  its  mission 
without  any  technical  difficulties.  This  limitation  can 
be  alleviated  by  introducing  breakdown  considerations 
similar  to  those  introduced  for  the  T.V. 

19.  Each  of  the  l  B.U.F.  requires  no  refuelling  during  the 
entire  operation. 


*  7 


20.  The  unloading  at  the  beach  may  be  sequential  or  in  para¬ 
llel,  but  not  both.  In  addition,  the  unloading  mode  at 
the  beach  is  independent  of  that  at  the  ship.  Moreover, 
in  the  parallel  unloading  mode  the  B.U.F.  performance  is 
not  influenced  by  the  number  of  B.U.F.  that  happens  to  be 
unloading  at  any  instant  of  time. 

21.  Times  TBC „  ~  .  depend  on  the  B.U.F.  characteristics  only, 
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and  are  independent  of  the  payload  unit  that  is  being  ser¬ 
viced  and  the  T.V.  from  which  the  payload  unit  in  question 
is  being  unloaded. 

22.  Time  TBC  8  ^  depends  on  the  B.U.F.  characteristics  only, 
and  is  independent  of  the  payload  unit  that  is  being  ser¬ 
viced. 

23.  The  payload  units  capable  of  self -induced  rolling  motion 
are  unloaded  from  the  T.V.  by  regular  B.U.F.  No  special 
provisions  have  been  made  in  this  program  to  include  any 
other  type  of  beach  unloading  facilities. 

24.  Each  of  the  L  B.U.F.  is  assumed  to  be  able  to  handle  the 
heaviest  and  bulkiest  payload  unit. 

25.  Only  one  payload  unit  may  be  handled  by  B.U.F.  at  any  one 
time. 

26.  Once  a  T.V.  beaches  at  the  Jlth  (l  =  1,2... L)  beach  unload¬ 
ing  area,  only  the  Jlth  B.U.F.  may  unload  it. 


27.  Each  of  the  I  T.V.  may  beach  at  one  beach  unloading  area 
only  during  the  execution  of  any  one  of  its  trips. 

28.  Each  of  the  L  B.U.F.  continues  to  unload  from  any  of  the 
I  T.V.  during  the  execution  of  any  one  of  its  trips  until 
all  the  payload  units  carried  by  the  T.V.  in  question  in 
that  particular  trip  have  been  unloaded. 

29.  When  the  ith  (i  =  li2...I)  T.V.  starts  from  its  base  it  is 
assumed  to  have  its  normal  fuel  tanks  fully  loaded  with 
fuel.  In  additional  tanks  that  are  so  constructed  as  to 
not  affect  the  volume  of  payload,  sufficient  fuel  is  car¬ 
ried  to  permit  the  ith  T.V.  to  reach  the  W.A.I  and  then 
the  appropriate  ship  unloading  area  where  it  hooks  up  and 
is  made  ready  for  the  unloading  operation,  without  using 
any  of  the  fuel  in  the  normal  fuel  tanks.  Furthermore,  it 
is  assumed  that  the  weight  of  this  additional  fuel  is 
smaller  than  or  equal  to  the  weight  of  the  payload. 

30.  The  fuel  stored  in  the  regular  fuel  tanks  must  be  suffi¬ 
cient  to  permit  any  T.V.  to  complete  at  least  one  round 
trip. 

31.  The  fuel  required  for  the  T.V.  to  return  to  its  base  after 
execution  of  Step  11,  as  described  in  the  mission  of  the 
ith  (i  =  2 1  2 ...  I )  T.V.  (see  Section  4)  is  less  than  or 
equal  to  the  fuel  required  for  it  to  complete  its  round 
trip. 


* 


32.  When  a  T.V.  is  to  be  refuelled,  it  is  always  refuelled 
completely . 

33.  The  refuelling  of  the  ith  (i  =  li2...I)  T.V.  is  to  be 
effected  only  when  it  is  alongside  the  mother  ship,  and  at 
whatever  ship  unloading  area  it  is  hooked  up. 

34.  The  fuel  requirements  of  the  T.V.  while  waiting  in  queues 
I  and  II,  while  hooked  up  alongside  the  mother  ship,  and 
while  beached  at  the  shore,  are  negligible. 

35.  The  time  T  .  depends  on  the  T.V.  characteristics  only, 
and  is  independent  of  the  ship  unloading  area  that  is  be¬ 
ing  serviced. 

36.  Each  of  the  I  T.V.  must  be  able  to  transport  the  heaviest 
and  bulkiest  payload  unit  in  any  of  its  trips . 

37.  The  T.V.  are  to  carry  payload  only  when  they  are  travel¬ 
ing  between  the  mother  ship  and  the  beach.  At  all  other 
times,  they  do  not  carry  any  payload. 

38.  The  mission  of  each  T.V.  is  such  that  it  does  not  inter¬ 
fere  with  that  of  any  other  T.V. 

39.  None  of  the  I  T.V.  requires  the  assistance  of  any  of  the 
S.U.F.  to  perform  Steps  3  and  6  as  described  in  the 
mission  of  the  ith  (i  =  1,2...  I)  T.V.  (see  Section  4). 


Sc 

40.  None  of  the  I  T.V.  requires  the  assistance  of  any  of  the 
B.U.F.  to  perform  Steps  9  and  11  as  described  in  the 
mission  of  the  ith  (i  =  li2...I)  T.V.  (see  Section  4). 

41.  The  times  T ^  1Q  jl-i  dePend  on  tiie  T*v*  characteristics 
only,  and  are  independent  of  the  S.U.F.  that  is  being 
serviced. 

42.  The  time  T Q  .  depends  on  the  T.V.  characteristics  only, 

v  j  t 

and  is  independent  of  the  S.U.F.  and  the  payload  units 
that  are  being  serviced. 

43.  The  times  _  .  depend  on  the  T.V.  characteristics 

only,  and  are  independent  of  the  B.U.F.  that  is  being 
serviced. 

44.  Each  T.V.  must  wait  (even  if  it  is  for  zero  time)  in 
queues  I  and  II  respectively  until  the  appropriate  unload¬ 
ing  areas  are  free  to  receive  them. 

45.  The  breakdown  considerations  for  processes  1,  3,  4,  6,  7 
and  8  (see  Table  3-1)  are  functions  of  i  only  and  are  time 
invariant. 

46.  The  breakdown  considerations  of  processes  2  and  5  (see 
Table  3-1)  are  functions  of  i  and  waiting  time  only  and 
are  again  time  invariant. 


47.  All  steps  (as  described  in  the  mission  of  the  ith 

(i  -•  li2...I)  T.V.)  that  can  be  performed  before  the 
breakdown  considerations  forcibly  remove  a  T.V.  from  our 
system,  are  executed.  (It  is  important  to  be  aware  of 
this  limitation  so  that  the  concept  of  antithetic  variance 
be  introduced  correctly  into  our  methodology.  It  can  be 
relaxed  at  will  by  changing  the  computer  program  accord¬ 
ingly.  ) 

48.  When  the  ith  (i  =  1,2...!)  T.V.  is  executing  any  other 
segment  of  its  trip  not  covered  by  processes  1-8  of 
Table  3-1,  it  is  assumed  that  it  cannot  malfunction. 

The  above  completes  our  discussion  of  the  fourth 
topic  of  the  Problem  Definition.  We  now  proceed  with  the  dis¬ 
cussion  of  the  final  topic  in  the  Problem  Definition,  namely, 
the  Selection  of  the  Figure  of  Merit. 


v)  Selection  of  the  Figure  of  Merit 


Finally,  as  was  already  stated  in  the  two  previous 
sections,  the  figure  of  merit  (the.  measure  of  success)  is  the 
weighted  combination  of  time  and  the  level  of  risk  and  un¬ 
certainty.  The  above  can  be  expressed  mathematically  as  fol¬ 
lows  : 


a  =  Minimize  (W  T  +  V)  (3.17/ 

where 

c  is  the  figure  of  merit, 

Wj  &  W 2  the  weighting  factors  such  that 
0  £  wyw2  £  - 
and  W i  +  W  g  =  1 , 

T  is  the  time  component  of  the  figure  of  merit 
involving  the  calculation  of  the  maximum  time 
elapsed  since  the  start  of  the  mission  to 

1.  prepare  the  Mother  Ship  for  departure  after 
all  the  payload  has  been  transferred  to  the 
T.V.  and  all  T.V.  have  cleared  the  Mother  Ship, 

2.  complete  the  payload  transfer  from  the  Mother 
Ship  to  point  A  on  the  beach, 

3.  return  all  T.V.  to  their  bases, 

4.  return  all  B.U.F.  to  their  bases, 

i.e.  T=Maximum  J TTM,AMAXSCc  ,  ,AMAXBC c  AMAXTV a  . I  (3.18) 

k=li  2.  .  .  K  V  6 8>1') 

Z—  1  j  2 .  •  .  L 
i=l j 2 ...  I 

U  is  the  level  of  risk  and  uncertainty  component  of 
the  figure  of  merit  involving  the  calculation  of 
the  percentage  of  T.V.  that  did  not  complete  their 
mission  (see  Section  4)  because  of  breakdown. 


and 


When  our  analysis  involves  no  breakdown  considera¬ 
tions,  then  U  in  equation  (3.17)  becomes  zero  and  our  figure 
of  merit  can  simply  be  written  as 

a  =  Minimize  (T) 

which  leads  directly  to  the  correct  solution. 


When  breakdown  considerations  are  involved  in  our 
analysis,  however,  equation  (3.17)  should  be  modified  in  or¬ 
der  to  give  the  correct  measure  of  success  for  each  case 
under  investigation.  This  modification  involves  the  normal¬ 
ization  of  the  decision  elements  T  and  U  to  T*  and  U'  ,  so 
that  the  normalized  decision  elements  have  the  same  effect  on 
our  decision. 


One  of  the  methods  that  can  accomplish  the  above  men¬ 
tioned  normalization  is  the  following: 


Let 


T  -  T  . 
_ mm 

T  -  T  . 
max  mm 


(3.19) 


and 


u1 


u  -  u, 


mm 


U  -  U  . 
max  mm 


(3.20) 


where 

and 


T  • 
mm 

U  . 
mm 


<  T  <  T 

<  U  <  U 


max 


max 


* 


(3.21) 

(3.22) 


where  Tmint  TmaxJ  Umin  and  Umax  can  be  determined  after  all 
the  cases  under  investigation  have  been  completed. 


♦When  no  breakdown  considerations  are  involved  in  our  analysis 


u  •  =  u  =  u 

mm  max 


0 


Thus  the  mathematical  form  of  the  figure  of  merit 
used  for  our  study  is 

a  -  Minimize  (W^  T  ’  +  W  ^  U')  (3.23) 

which  always  leads  to  the  desired  solution  regardless  of 
whether  breakdown  considerations  are  included  in  our  study 
or  not. 

The  above  completes  the  discussion  on  the  Problem 
Definition,  and  now  we  may  proceed  with  the  discussion  on  the 
formulation  of  the  Mathematical  Model  for  our  system. 


vi)  Breakdown  Considerations 


If  IDBRTV  .  .  =  1  (ja  =  1,2...  8;  i  -  1,2... I)  there  are 

no  breakdown  considerations  involved  in  our  analysis  for  the 
part  of  the  journey  of  the  ith  T.V.  described  by  the  j'^th 
process,  and  so  the  discussion  given  here  does  not  apply  for 
that  T.V.  and  that  part  of  its  journey. 


Before  the  ith  T.V.  is  allowed  to  complete  its  journey  to 
W.A.I  from  its  base,  a  random  number,  R,  is  generated  from  the 
appropriate  distribution,  as  dictated  by  IDBRTV  .,  utilizing 

i,i 

INBRTV .  .  as  the  first  seed*.  If 

i  j  “t 


R  <  BRKTV 


l,i 


(4.1) 


then  the  ith  T.V.  is  assumed  lost  from  our  system.  If 
equation  (4.1)  is  not  satisfied,  then  the  ith  T.V.  is  allowed 
to  enter  W.A.I  and  it  is  assumed  that  no  breakdown  has 
occurred  in  the  part  of  ,tho  journey  of  the  ith  T.V.  described 
above . 


In  the  mathematical  model  of  our  system,  a  record  of  the 
T.V.  lost  at  this  stage  is  kept.  If,  at  the  end  of  the  Mother 
Ship's  mission,  any  of  these  T.V.  has  not  been  selected  at  all 
by  the  appropriate  T.V.  use  strategy  for  W.A.I,  a  note. to  the 
user  is  given  stating  that  fact. 


Before  the  ith  T.V.  is  allowed  to  depart  from  W.A.I,  a 
random  number,  R,  is  generated  from  the  appropriate 

*The  updated  INBRTV .  .are  used  as  the  seeds  of  subsequent 
generations .  * v 


distribution,  as  dictated  by  IDBRTV  0  .  utilizing  INBRTV .  . 

"  t  U  j  I, 

as  the  first  seed*. 

If 

TTVP.  . 

R  <  - LL L  •*  BRKTV-  .  (4.2) 

WT1MAX  ^ 

where  TTVP  6  ^  =  total  waiting  time  of  the  ith  T.V.  in  W.A.I 
J  in  this  trip  of  the  T.V.  in  question. 

if  TTVP-  .  <  WT1MAX  (4.3) 

0  j  2* 

or  if  R  <  BRKTV -  .  (4.2) 

2  s  X 

if  TTVP  .  >  WT1MAX  (4.3) 


then  the  ith  T.V.  is  assumed  lost  from  our  system.  If  equa¬ 
tion  (4.2)  is  not  satisfied,  then  the  ith  T.V.  is  allowed  to 
depart  from  W.A.I  for  the  appropriate  ship  unloading  area  and 
it  is  assumed  that  no  breakdown  has  occurred  while  the  T.V.  in 
question  was  waiting  j”  W.A.I. 


Before  the  ith  T.».  is  allowed  to  reach  the  appropriate 
ship  unloading  area  from  W.A.I,  a  random  number,  R,  is  generat¬ 
ed  from  the  appropriate  distribution,  as  dictated  by  IDBRTV _ 

Vj  t 

utilizing  INBRTV s  ^  as  the  first  seed**. 

If 

R  <  BRKTV _  .  (4.4) 

O  j  2> 

then  the  ith  T.V.  is  assumed  lost  from  our  system.  If  equation 
(4.4)  is  not  satisfied  then  the  ith  T.V.  is  allowed  to  reach 


*The  updated  INBRTV 2  i  are  used  as  the  seeds  of  subsequent  gen 
erations.  * 

**The  updated  INBRTV _  .  are  used  as  the  seeds  of  subsequent 
generations .  * 1 


the  ship  unloading  area  in  question  and  it  is  assumed  that 
no  breakdown  has  occurred  in  the  part  of  the  trip  of  the  ith 
T.V.  described  above. 


Before  the  ith  T.V.  is  allowed  to  reach  W.A.II  from  the 
appropriate  ship  unloading  area,  a  random  number ,  R,  is  gene¬ 
rated  from  the  appropriate  distribution,  as  dictated  by 
IDBRTV ,  ,t  utilizing  INBRTV .  .  as  the  first  seed*. 

*  j  *  *  J  «' 

If  R  <  BRKT .  .  (4.5) 

then  the  ith  T.V-  is  assumed  lost  from  our  system.  If  equa¬ 
tion  (4.5)  is  not  satisfied,  then  the  ith  T.V.  is  allowed  to 
enter  W.A.II  and  it  is  assumed  that  no  breakdown  has  occurred 
in  the  part  of  the  trip  of  the  ith  T.V.  described  above. 


Before  the  ith  T.V.  is  allowed  to  depart  from  W.A.II,  a 
random  number,  R,  is  generated  from  the  appropriate  distribu¬ 
tion,  as  dictated  by  IDBRTV,.  .,  utilizing  INBRTV c  .  as  the 

*  t>j  %  t>  j  U 

first  seed**. 

TTVP.  . 

If  R  <  *  BRKTV  .  (4.6) 

WT2MAX 

where  TTVP„  .  *  total  waiting  time  of  the  ith  T.V.  in 
J  W.A.II  in  this  .trip  of  the  ith  T.V. 

if  TTVP.  .  <  WT2MAX  (4.7) 

Oj  t 

or  if  R  <  BRKTV e  .  (4.6) 

if  TTVP.  .  >  WT2MAX  (4.7) 

o,  t  — 

♦The  updated ' INERT V~.  .  are  used  as  the  seeds  of  subsequent 
generations .  * 1 

**The  updated  INBRTV .  .  are  used  as  the  seeds  of  subsequent 
generations .  J  % 


then  the  ith  T.V.  is  assumed  lost  from  our  system.  If 
equation  (4.6)  is  not  satisfied,  then  the  ith  T.V.  is  allowed 
to  depart  from  W.A.II  for  the  appropriate  beach  unloading  area 
and  it  is  assumed  that  no  breakdown  has  occurred  while  the 
T.V.  in  question  was  waiting  in  W.A.II. 

Before  the  ith  T.V.  is  allowed  to  reach  the  appropriate 
beach  unloading  area  from  W.A.II,  a  random  number,  R,  is 
generated  from  the  appropriate  distribution,  as  dictated  by 
IDBRTV „  .,  utilizing  IRBRTV *  .  as  tie  first  seed*.  If 

o  ,  1  b  j  t 

R  <  BRKTV.  .  (4.8) 

o  ,  i 

then  the  ith  T.V.  is  assumed  lost  from  our  system.  If  eTua- 
tion  (4.8)  rs  not  satisfied,  then  the  ith  T.V.  is  allowed  to 
reach  the  beach  unloading  area  in  question  and  it  is  assumed 
that  no  breakdown  has  occurred  in  the  part  of  the  trip  of  the 
ith  T.V.  described  above. 

Before  the  ith  T.V.  is  allowed  to  reach  W.A.I  from  the 
appropriate  beach  unloading  area,  a  random  number,  R,  is  gen¬ 
erated  from  the  appropriate  distribution,  as  dictated  by 
IDBRTV 7  utilizing  IRBRTV 7  .  as  the  first  seed**.  If 

R  <  BRKTV7  .  (4.9) 

*  The  updated  IRBRTV *  .  are  used  as  the  seeds  of  subsequent 
generations .  * 

**The  updated  IRBRTV 7  .  are  used  as  the  seeds  of  subsequent 
generations .  * 


then  the  £th  T,V*  is  assumed  lost  from  our  system.  If 
equation  (4.9)  is  not  satisfied,  then  the  £th  T.V.  is  allowed 
to  reach  W.A.I  and  it  is  assumed  that  no  breakdown  has  occurred 
in  the  part  of  the  trip  cf  the  i th  T.V.  described  above. 

In  our  mathematical  model  the  action  taken  as  described 
above  is  forfeited  for  all  the  T.V.  that  happen  to  be  i.i. 

W.A.I  immediately  after  the  appropriate  T.V.  departed  from 
W.A.I  to  load  the  last  payload  unit(s)  from  the  Mother  Ship. 

All  these  T.V.  are  assumed  to  have  started  for  their  res^-ect- 
ive  bases  at  the  times  when  they  last  started  the  execution  of 
Step  11  given  in  the  mission  of  the  ith  (i  =  lt2..,I)  T.V. 


Before  the  £th  T.V.  is  allowed  to  reach  its  base  from  the 
appropriate  beach  unloading  area,  a  random  number,  R,  is  gen¬ 
erated  from  the  appropriate  distribution,  as  dictated  by 
IDBRTV 0  .,  utilizing  INBRTV  0  .  as  the  first  seed*.  If 

O  jV  o j V 

R  <  BRKTV 0  .  (4.10) 

o  ,  v 

then  the  £th  T.V.  is  assumed  lost  from  our  system.  If  equa¬ 
tion  (4.10)  is  not  satisfied,  then  the  ith  T.V.  is  allowed  to 
reach  its  base  and  it  is  assumed  that  no  breakdown  has  oc¬ 
curred  in  the  part  of  the  journey  of  the  ith  T.V.  described 
above. 


* 


The  updated  INBRTVg  .  are  used  as  th.>  seeds  of  subsequent 
generations.  J 


The  payload  units  are  first  separated  into  K  groups  of 
(k  =  lt2...K )  units  .each.  Next,  the  payload  units  of  the 
first  group  are  identified  by  distinct  and  ascending  numbers. 
The  numbers  start  from  1  and  finish  at  N  .  Then  the  payload 
units  of  the  second  group  are  similarly  identified  by  numbers 
starting  fromf  1+N  2)  to  In  the  same  manner,  the  pay¬ 

load  units  of  the  kth.  group  are  identified  by  numbers  from 
k-1  k  1 

1+  l  N  .  to  l  N  .  and  so  on,  until  the  last  payload  unit  of 
{  3-1  3 ]  (3-1  3 j  [{  ’ 

the  last  group  is  identified  as  the  l  N  •  th  payload  unit. 

U-i 3  J 

Once  the  above  mentioned  identification  is  complete,  the  pay- 

load  units  are  loaded  into  the  Mother  Ship  in  the  following 

manner.  For  each  group*,  the  payload  unit  with  the  smallest 

identification  number  is  loaded  first  into  the  appropriate 

location  in  the  Mother  Ship.  After  that,  the  th  payload 

unit  of  each  of  the  K  groups  is  loaded  immediately  after  the 

payload  unit  of  the  same  group  into  the  appropriate 

k-1  k 

location  in  the  Mother  Ship,  where  2  +  I  N  •  i  J*.  <  l  N .  . 

0=1  0  o-l  3 

This  operation  continues  until  all  payload  units  have  been 
loaded  into  the  Mother  Ship. 


*Note  that  the  loading  operations  of  each  of  the  K  groups  are 
totally  independent  of  one  another.  The  only  point  of 
interest  to  us  is  the  loading  order  of  the  payload  units  of 
each  of  the  K  groups. 


The  unloading  strategy  we  adopted  for  this  study  is 
the  following.  The  payload  units  of  each  of  the  K  groups* 
are  to  be  unloaded  in  the  exact  reverse  of  the  order  in  which 
they  were  loaded  into  the  Mother  Ship.  Also,  the  payload 
units  are  to  be  unloaded  at  the  beach  in  the  exact  reverse  of 
the  order  in  which  they  were  loaded  in  the  T.V.  in  question. 


*The  unloading  sequences  of  each  of  the  K  groups  are  totally 
independent  of  one  another. 


a  i 
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viii)  l.V.  Use  Strategy  A  for  W.A.I. 

This  use  strategy  (see  subroutine  ASLTVA)  is  the  first  of 

the  strategies  incorporated  in  the  mathematical  model  for 

selecting  T.V.  from  W.A.I.  It  is  basically  a  first  come,  first 

served  strategy  and  is  oriented  tovzards  volume- limited  T.V. 

The  T.V.  selection  is  governed  by  the  rollowing  rules: 

1.  First  come,  first  selected*.  In  the  event  of  a  tie 
select  the 

2.  T.V.  with  the  biggest  (available  voljune/time) ** .  In  the 
event  of  a  further  tie  select  the 

3.  T.V.  with  the  biggest  (available  weight  capacity/ time) ***. 

In  the  event  of  a  further  tie  select  the 

4.  T.V.  with  the  biggest  available  volume.  In  the  event  of  a 
farther  tie  select  the 

5.  T.V.  with  the  biggest  available  weight  capacity.  In  the 
event  of  a  further  tie  select  the 

6.  Speedier****  T.V.  In  the  event  of  a  further  tie  select  the 

7.  T.V.  with  the  smallest  identification. 


*Not  necessarily  in  effect  if  the  T.V.  selected  cannot  be 
serviced  immediately.  All  T.V.  available  in  W.A.I  before  the 
appropriate  ship  unloading  area  is  free  are  equally  eligible 
as  far  as  the  first  rule  is  concerned. 

♦♦Available  volume  specified  by  AMAXTV 7  ±  (i=l,2. . .1) ,  time  as 
measured  by  ( TTV  4 ^  +  TTV2o,i  +  TTV22/i2  i=l,2. .  .1) . 

♦♦♦Available  weight  capacity  specified  by  AMAXTV .  (1=1,2. .  .1) 

****Speed  measured  by  the  above  mentioned  time,  noting  that  the 
smaller  (TTV4  v  +  TTV 10 ^  +  TTVn  a  i=l,2...I )3  then  the 
speedier  the  T;V.  is. 


Co  3 


ix)  T. V.  Use  Strategy  B  for  W.A.I. 

This  use  strategy  (see  subroutine  ASLTVB )  is  the  second 
of  the  strategies  incorporated  in  the  mathematical  model  for 
selecting  T.V.  from  W.A.I. 

It  is  again  basically  a  first  come/  first  served  strategy 
but  is  oriented  towards  weight-limited  T.V. 

The  T.V.  selection  is  governed  by  the  same  rules  given 
for  the  T.V.  use  strategy  A  for  W.A.I,  with  the  exception  that 
the  second  rule  is  inter changed  with  the  third,  and  the  fourth 
with  the  fifth. 


x)  T.V.  Use  Strategy  A  for  W.A.II. 

This  use  strategy  (see  subroutine  BSLTVA)  is  the  first 
of  the  strategies  incorporated  into  the  mathematical  model  for 
selecting  T.V.  from  W.A.II*.  It  is  basically  a  first  come, 
first  served  strategy  and  is  oriented  towards  volume-limited 
T.V.  (For  this  reason,  it  is  logical  that  this  use  strategy 
should  be  employed  for  the  T.V.  selection  from  W.A.II  when  use 
strategy  A  for  W.A.I  is  employed  to  select  T.V.  from  W.A.I.) 

The  T.V.  selection  is  governed  by  the  following  rules: 

1.  First  come,  first  selected**.  In  the  event  of  a  tie  select 
the 

2.  T.V.  with  the  biggest  (available  volume/time) *** .  In  the 
event  of  a  further  tie  select  the 

3.  T.V.  with  the  biggest  (available  weight  capacity /time) **** . 
In  the  event  of  a  further  tie  select  the 

4.  T.V.  with  the  biggest  available  volume.  In  the  event  of  a 
further  tie  select  the 

5.  T.V.  with  the  biggest  available  weight  capacity.  In  the 
event  of  a  further  tie  select  the 

*The  T.V.  selected  from  W.A.I  is  retained  until  a  T.V.  is 
available  at  an  earlier  time  in  W.A.II. 

**Not  necessarily  in  effect  if  the  T.V.  selected  cannot  be 
serviced  immediately.  All  T.V.  available  in  W.A.II  before  the 
appropriate  beach  unloading  area  is  free  are  equally  eligible 
as  far  as  the  first  rule  is  concerned. 

***Available  volume  specified  by  AMAXTV?  ^  (i=l,2. . . I) ,  time 
as  measured  by  (TTV1S* ^  +  TTV +  TTV\6^^s  i=l,2. . .1) . 

****Available  weight  capacity  specified  by  AMAXTVg  i 

(1=1,2...!). 


6.  T.V.  with  the  biggest  (utilized  volume/time}* .  In  the 
event  of  a  further  tie  select  the 

7.  T.V.  with  the  biggest  (utilized  weight  capacity/time) **. 

In  the  event  of  a  further  tie  select  the 

8.  T.V.  with  the  biggest  utilized  volume.  In  the  event  of  a 
further  tie  select  the 

9.  T.V.  with  the  biggest  utilized  weight  capacity.  In  the 
event  of  a  further  tie  select  the 

10.  Speedier***  T.V.  In  the  event  of  a  further  tie  select  the 

11.  T.V.  with  the  smallest  identification. 


♦Utilized  volume  is  measured  by  summing  up  the  volume  of  each 
payload  unit  and  its  lashings  carried  by  the  T.V.  in  question 
in  that  particular  trip,  time  as  measured  by  (TTV1,  .  + 

ITV.C  .  +  TTV,a  i=lt 2. . . I) .  ‘ 

1  O  j  'V  i  0  j 

♦♦Utilized  weight  capacity  is  measured  by  summing  up  the 
weight  of  each  payload  unit  and  its  lashings  carried  by  the 
T.V.  in  question  in  that  particular  trip. 


♦♦♦Speed  measured  by  the  abo’” 
smaller  (TTV,,  .  +  TTV n ,  . 


antioned  time,  noting  that  the 

.),  then  the  speedier  the 

1  0  •  *1* 
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xi )  T.V.  Use  Strategy  B  for  W.A.II. 


This  use  strategy  (see  subroutine  BSLTVB )  is  the  second 
of  the  strategies  incorporated  into  the  mathematical  model 
for  selecting  T.V.  from  W.A.II*. 

It  is  again  basically  a  first  come,  first  served  strategy 
but  is  oriented  towards  weight-limited  T.V.  (For  this  reason, 
it  is  logical  that  this  use  strategy  should  be  employed  for 
the  T.V.  selection  from  W.A.II  when  use  strategy  B  for  W.A.I 
is  employed  to  select  T.V.  from  W.A.I.) 

The  T.V.  selection  is  governed  by  the  same  rules  given 
for  the  T.V,  use  strategy  A  for  W.A.II  with  the  exception  that 
the  second  rule  is  interchanged  with  the  third,  the  fourth 
with  the  fifth,  the  sixth  with  the  seventh,  and  the  eighth 
with  the  ninth. 


♦The  T.V.  selected  from  W.A.I  is  retained  until  a  T.V.  is 
available  at  an  earlier  time  in  W.A.II. 


.f 


xii)  S.U.F.  Use  Strategy  A 

This  use  strategy  (see  subroutine  SLSCA )  is  the  first  of 

the  strategies  incorporated  into  the  mathematical  model  for 

selecting  S.U.F.  It  is  a  first  come,  first  served  strategy. 

The  S.U.F.  selection  is  governed  by  the  following  rules: 

1.  The  S.U.F.  whose  associated  ship  unloading  area  is  first 
free  is  first  selected.  In  the  event  of  a  tie,  select  the 

2.  S.U.F.  that  has  the  most  cargo  to  unload  at  the  time  the 
above  mentioned  tie  occurred.  In  the  event  of  a  further 
tie  select  the 

3.  S.U.F.  that  will  be  ready  to  commence  its  unloading  cycle 
first*.  In  the  event  of  a  further  tie  select  the 

4.  Speedier**  S.U.F.  In  the  event  of  a  further  tie  select  the 

5.  S.U.F.  with  the  smallest  identification. 


*A.  S.U.F.  is  ready  to  commence  its  unloading  cycle  at  the 
instant  it  reaches  the  ship  unloading  area  associated  with  it 
for  the  first  time  or  at  any  other  time  immediately  after  the 

S. U.F,  in  question  is  made  ready  to  travel  again  after  it  has 
unloaded  the  appropriate  payload  unit  into  the  appropriate 

T. V. 

**Speed  measured  by  the  time  (TSCz^  +  TSC3  ^  +  TSC4  £  + 

+  TSC?^;  k=lt2. . .K) ,  noting  that  tAe  smaller  this 
time^  then  the  speedier  the  S.U.F.  is. 


* 
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xiii)  S.U.F.  Use  Strategy  B . 

This  use  strategy  (see  subroutine  SLSCB)  is  the  second 
of  the  strategies  incorporated  into  the  mathematical  model 
for  selecting  S.U.F.  It  is  basically  a  first  come,  first 
served  strategy. 

The  S.U.F.  selection  is  governed  by  the  same  rules  given 
for  S.U.F.  use  strategy  A,  with  the  exception  that  Rule  1  is 
not  necessarily  in  effect  if  the  S.U.F.  has  to  await  the 
arrival  of  a  T.V.  in  W.A.I.  In  that  case,  all  S.U.F.  whose 
associated  ship  unloading  areas  are  free  before  the  arrival  of 
the  T.V.  in  question  in  W.A.I  are  equally  eligible  regarding 


the  first  rule. 


This  use  strategy  (see  subroutine  SLBCA)  is  the  first  of 


the  strategies  incorporated  into  the  mathematical  model  for 
selecting  B.U.F.  It  is  a  first  come,  first  served  strategy. 
The  B.U.F.  selection  is  governed  by  the  following  rules: 

1.  The  B.U.F.  whose  associated  beach  unloading  area  is  first 
free  is  first  selected.  In  the  event  of  a  tie  select  the 

2.  B.U.F.  that  will  be  ready  to  commence  its  unloading  cycle 
first*.  In. the  event  of  a  further  tie  select  the 

3.  Speedier**  B.U.F.  In  the  event  of  a  further  tie  select  the 

4.  B.U.F.  with  the "smallest  identification. 


*B.U.F.  is  ready  to  commence  its  unloading  cycle  at  the  in¬ 
stant  it  is  made  ready  for  travel  after  it  has  reached  point 
A  on  the  beach  for  the  first  time  or,  at  any  other  time, 
after  the  appropriate  payload  unit  has  been  released  at  point 
A  on  the  beach  and  the  B.U.F.  has  been  made  ready  to  travel 
again. 

**Speed  measured  by  the  time  (TBC4  ^  +  TBC6  ^ +TBC?  Si+ TBCg 
1*1,2. .  .L)  noting  that  the  smaller  the  albove  mentioned  time, 
then  the  speedier  the  B.U.F.  is. 


This  use  strategy  (see  subroutine  SLBCB )  is  the  second  of 
the  strategies  incorporated  into  the  mathematical  model  for 
selecting  B.U.F.  It  is  basically  a  first  come,  first  served 
strategy. 

The  B.U.F.  selection  is  governed  by  the  same  rules  given 
for  B.U.F.  use  strategy  A,  with  the  exception  that  Rule  1  is 
not  necessarily  in  effect  if  the  B.U.F.  has  to  await  the 
arrival  of  a  T.V.  in  W.A.II.  In  that  case,  all  B.U.F.  whose 
associated  beach  unloading  areas  are  free  before  the  arrival 
of  the  T.V.  in  question  in  W.A.II  are  equally  eligible  regard¬ 
ing  the  first  rule. 


The  above  completes  our  description  of  the  formulation  of 
the  mathematical  model  for  our  study.  We  now  proceed  with  a 
discussion  of  the  solution  method. 


4.  Formulation  of  the  Mathematical  Model 

Instead  of  presenting  the  mathematical  equations  that 
make  up  the  mathematical  model  for  our  study,  the  mission  of 
each  subsystem  involved,  the  breakdown  considerations ,  and  the 
description  of  the  use  strategies  are  given.  The  reason  for 
choosing  this  approach  in  describing  the  formulation  of  our 
mathematical  model  is  because  we  achieve  our  goal  more  easily, 
as  the  nature  of  the  underlying  mathematical  equations  makes 
them  very  cumbersome  to  write.  If  the  user  desires  to  write 
the  mathematical  equations  of  the  model,  he  can  do  so  by  refer 
ring  to  the  computer  code  appended  at  the  end  of  this  study, 
or  by  translating  the  descriptions  given  here  in  the  form  of 
mathematical  equations. 

As  was  stated  in  the  second  section  of  this  report,  the 
mathematical  model  of  our  study  is  so  formulated  as  to  be  an 
exact  replica  of  the  system  under  consideration  to  the  degree 
of  accuracy  and  extent  desired.  Also,  because  of  the  type  of 
problem  we  are  dealing  with,  the  nature  of  the  model  is 
stochastic.  In  setting  up  the  mathematical  model,  care  was 
taken  to  keep  it  as  simple  as  possible  to  permit  easy  analysis 
and  yet  to  construct  it  so  that  it  exhibits  all  the  phenomena 
under  consideration,  as  required. 

The  description  of  the  formulation  of  our  mathematical 
model  is  presented  in  the  following  manner: 


i)  Description  of  the  mission  of  the  Mother  Ship. 

ii)  Description  of  the  unloading  procedure  of  the  nth 

K 

(n  =  1,2...  I  N.)  payload  unit. 
k=l  K 

iii)  Description  of  the  mission  of  the  kth  (k  *  1,2...K) 
S.U.F. . 

iv)  Description  of  the  mission  of  the  Jtth  (Si  =  1,2...L) 
B.U.F. 

v)  Description  of  the  mission  of  the  ith  (i  =  1,2...I)  T. 

vi)  Description  of  the  breakdown  considerations, 

vii)  Description  of  the  payload  unloading  strategies, 

viii)  Description  of  the  T.V.  use  strategy  A  for  W.A.I. 

ix)  Description  of  the  T.V.  use  strategy  B  for  W.A.I. 

x)  Description  of  the  T.V.  use  strategy  A  for  W.A.II. 

xi)  Description  of  the  T.V.  use  strategy  B  for  W.A.II. 

xii)  Description  of  the  S.U.F.  use  strategy  A. 

xiii)  Description  of  the  S.U.F.  use  strategy  B. 

xiv)  Description  of  the  B.U.F.  use  strategy  A. 

xv)  Dsscription  of  the  B.U.F.  use  strategy  B. 
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i)  Mission  of  the  Mother  Ship. 

Step  1.  The  Mother  Ship  arrives  at  the  theater  of  operations. 

Step  2.  The  Mother  Ship  is  properly  moored  in  position  after 
it  arrives  at  the  theater  of  operations.  Upon  com¬ 
pletion  of  this  step,  all  ship  unloading  areas  become 
free  for  the  first  time. 

Step  3.  The  Mother  Ship  is  freed  from  its  moorings  and  is 

made  ready  to  travel  after  all  S.U.F.  are  secured  in 
position  and  all  T.V.  have  cleared  the  Mother  Ship. 

Upon  completion  of  Step  3  given  above,  the  Mother  Ship's 

mission  is  completed. 


K 

ii)  Unloading  Procedure  of  the  nth  (n  =  1,2...  I  N.) 

Payload  Unit. _ k=l 

Step  1  The  nth  payload  unit  is  released  after  the  appropriate 
S.U.F.  has  reached  the  payload  unit  in  question. 

Step  2  The  nth  payload  unit  is  secured  on  the  appropriate 
S.U.F.  after  the  payload  unit  in  question  has  been 

\ 

released. 

Step  3  The  nth  payload  unit  is  transported  to  the  appropriate 
ship  unloading  area  after  it  has  been  secured  on  the 

S. U.F.  in  question. 

Step  4  The  nth  payload  unit  is  unloaded  into  the  appropriate 

T. V.  and  then  freed  from  the  S.U.F.  in  question.  This 
operation  is  performed  only  if 

a)  the  appropriate  T.V.  has  been  properly  secured  in 
the  ship  unloading  area  and  has  been  made  ready 
for  the  loading  operation  and  has  completed  its 
refueling  (if  refueling  was  necessary) , 

b)  the  previous  payload  unit  unloaded  by  the  S.U.F. 
in  question  is  fully  secured  in  the  T.V.  in. ques¬ 
tion  (this  requirement  is  void  if  the  nth  payload 
unit  is  the  first  payload  unit  to  be  unloaded  in 
any  of  the  T.V. 's  trips),  and 

c)  the  T.V. 's  remaining  capacity  can  accept  the  nth 


payload  unit. 


If  that  is  not  the  case  the  nth  payload  unit  will 
have  to  wait  in  the  appropriate  ship  unloading  area 
until  the  above  requirements  are  satisfied. 


Step  5  The  nth  payload  unit  is  secured  in  position  on  the 
T.V.  in  question  after  it  has  been  freed  from  the 
appropriate  S.U.F. 

Step  6*  The  nth  payload  unit  is  transported  to  W.A.II  from  the 
appropriate  ship's  unloading  area.  This  operation  is 
performed  only  if 

a)  the  capacity  of  the  T.V.  in  question  is  such  that 
it  cannot  accept  the  next  payload  unit  or 

b)  the  nth  payload  unit  was  the  last  payload  unit  as¬ 
sociated  with  the  S.U.F.  in  question.  If  that is  not 
the  case,  the  nth  payload  unit  will  wait  alongside 
the  ship  unloading  area  in  question  until  one  or 
both  of  the  above  requirements  are  satisfied.  In 
addition,  please  note  that  the  operation  involved 
in  Setp  6  is  performed  after 

a)  the  last  payload  unit  unloaded  into  the  T.V.  in 
question  in  any  of  its  trips  has  been  properly 
secured,  and 

b)  after  the  T.V.  in  question  has  been  unhooked 
from  the  appropriate  ship  unloading  area  and  has 
been  made  ready  for  travel. 

♦Steps  6-ll  will  not  be  executed  if  the  T.V.  in  question  breaks 
down  while  on  route  to  W.A.II  from  the  Mother  Ship.  The  nth 
payload  unit  together  with  the  T.V.  in  question  is  then 
assumed  lost  from  our  system. 


Step  7*  The  nth  payload  unit  waits  in  W.A.II  after  it  has 

arrived  there  until  the  appropriate  B.U.F.  is  ready 
to  receive  the  T.V.  that  carries  the  payload  unit 
in  question. 

Step  8**The  nth  payload  unit  is  transported  to  the  appro¬ 
priate  beach  unloading  area  from  W.A.II  after  it  has 
waited  there  appropriately  until  '.he  beach  unloading 
area  in  question  is  freed. 

Step  9  The  nth  payload  unit  is  released,  utilizing  the 

appropriate  B.r.F.  The  operation  is  performed  only  if 
the  T.V.  in  question  is  beached  and  ready  to  commence 
the  unloading  operation,  and  the  appropriate  B.U.F. 
has  reached  the  appropriate  beach  unloading  area.  If 
that  is  not  the  case,  the  releasing  of  the  nth  payload 
unit  has  to  wait  until  the  above  requirements  are 
satisfied. 

Step  10  The  nth  payload  unit  is  secured  on  the  appropriate 
B.U.F.  after  it  has  been  released. 

1  ■ 

Steps  7-11  will  not  be  executed  if  the  T.V.  in  question  breaks 
down  while  waiting  in  W.A.II.  The  nth  payload  unit  together 
with  the  T.V.  in  question  is  then  assumed  lost  from  our 
system. 

**Steps  8-11  will  not  be  executed  if  the  T.V.  in  question 
breaks  down  while  on  route  to  the  beach  from  W.A.II.  The 
nth  payload  unit  together  with  the  T.V.  in  question  is  then 
assumed  lost  from  our  system. 


.* 
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Step  11  After  the  nth  payload  unit  is  secured  on  the  appro¬ 
priate  B.U.F.  it  is 

a)  transported  from  the  appropriate  beach  un¬ 
loading  area  to  point  A  on  the  beach  by 
utilizing  the  appropriate  B.U.F. , 

b)  unloaded  at  point  A,  and  then 

c)  freed  from  the  B.U.F.  in  question. 


Upon  completion  of  Step  11  given  above,  the  nth  payload 


unit  is  considered  to  have  reached  its  destination. 


iii)  Mission  of  the  kth  (k 


2...K)  S.U.F 


Step  1  The  kth  S.U.F.  is  made  ready  to  start  the  unloading 
operation  and  is  allowed  to  reach  the  kth  ship  un¬ 
loading  area  after  the  Mother  Ship  is  properly 
moored. 

Step  2  The  kth  S.U.F.  travels  to  the  appropriate  payload  unit 
from  the  kth  ship  unloading  area  after  it  has  been 
made  ready  to  travel,  or  after  the  above  step  is  com¬ 
pleted. 

Step  3  The  kth  S .U.F>^eleases  the  appropriate  payload  unit 
after  the  kth  S.U.F>^ias  reached  the  payload  unit  in 
question.  \ 

Step  4  The  kth  S.U.F.  has  the  payload  unit  in  question 

secured  onto  it  after  the  payload  unit  in  question  has 
been  released. 

The  kth  S.U.F.  transports  the  payload  unit  in  question 
to  the  kth  ship  unloading  area  after  the  above  men¬ 
tioned  payload  unit  has  been  secured  on  the  kth  S.U.F. 

\ 

The  kth  S.U.F.  unloads  the  payload  unit  in  question 
into  the  appropriate  T.V.  and  releases  it,  and  then 
the  kth  S.U.F.  is  made  ready  to  travel  again.  This 
operation  is  performed  only  if 


a)  the  appropriate  T.V.  is  properly  secured  in  the 
fcth  ship  unloading  area  and  has  been  made  ready  for 
the  loading  operation  and  has  completed  its  refuel¬ 
ing  (if  refueling  was  necessary) , 

b)  the  previous  payload  unit  unloaded  by  the  fcth 
S.U.F.  is  fully  secured  in  the  T.V.  in  question 
(this  requirement  is  void  if  the  payload  unit  in 
question  is  the  first  payload  unit  to  be  unloaded 
in  any  of  the  T.V. 's  trips),  and 

c)  the  T.V. 's  remaining  capacity  can  accept  the  pay- 
load  unit  in  question. 

If  any  of  the  above  is  not  satisfied,  the  fcth  S.U.F. 
must  wait  until  all  three  requirements  given  above  are 
satisfied. 

Step  7  •  The  fcth  S.U.F.  is  allowed  to  travel  back  to  its  origin¬ 
al  position  from  the  fcth  ship  unloading  area  and  then 
be  secured  to  its  original  position  after  all  the 
payload  units  have  been  unloaded  into  the  T.V. 

Upon  completion  of  Step  7  given  above,  the  mission  of  the 
fcth  S.U.F.  is  completed.  Please  note  that  when  Steps  2-6  given 
above  are  executed,  the  fcth  S.U.F.  is  then  said  to  have  execut¬ 
ed  one  unloading ■ cycle.  It  then  follows  that  during  its  mis¬ 
sion,  the  fcth  S.U.F.  will  execute  unloading  cycles.  The 
completion  of  an  unloading  cycle  automatically  starts  the  next 
one  until  all  ^  unloading  cycles  have  been  executed. 


iv)  Mission  of  the  Jtth  CSL  =  B.U.F. 

Step  1  The  Jtth  B.U.F.  starts  from  its  base  having  as  its  fin¬ 
al  destination  point  A  on  the  beach. 

Step  2  The  Jtth  B.U.F.  is  made  ready  to  start  the  unloading  op¬ 
eration  after  the  Jtth  B.U.F.  has  reached  point  A  on 
the  beach.  Upon  completion  of  this  step/  all  beach 
unloading  areas  become  free  for  the  first  time. 

Step  3  The  Jtth  B.U.F.  travels  to  the  Jtth  beach  unloading  area 
after  it  has  been  made  ready  to  travel. 

Step  4  The  Jtth  B.U.F.  releases  the  appropriate  payload  unit 
after  it  has  reached  the  Jtth  beach  unloading  area. 

This  operation,  of  course,  is  performed  only  if  the 
appropriate  T.V.  is  beached  and  is  ready  to  commence 
the  unloading  operation.  If  that  is  not  the  case,  the 
Jtth  B.U.F.  must  wait  for  the  T.V.  in  question  to  ar¬ 
rive,  to  be  beached  and  to  be  made  ready  for  the  un¬ 
loading  operation. 

Step  5  The  Jtth  B.U.F.  has  the  payload  unit  in  question  se¬ 
cured  onto  it  after  the  payload  unit  has  been  released. 

Step  6  The  Jtth  B.U.F.  (after  the  payload  unit  in  question  has 
been  secured  onto  it)  transports  the  payload  unit  in 
question  from  the  Jtth  beach  unloading  area  to  point  A 
on  the  beach.  After  that,  it  releases  the  unit  and  the 


B.U.F.  in  question  is  made  ready  to  travel  again. 


Step  7  The  Ath  B.U.F.  is  secured  in  its  original  position  and 
is  prepared  for  departure  after  its  unloading  mission* 
has  been  completed. 

Step  8  The  Ath  B.U.F.  is  allowed  to  return  to  its  base  after 
it  has  been  prepared  for  departure. 

Upon  completion  of  Step  8  given  above,  the  mission  of  the 
Ath  B.U.F.  is  completed.  Please  note  that  when  Steps  3-6 
given  above  are  executed,  the  Ath  B.U.F.  is  said  to  have  exe¬ 
cuted  one  unloading  cycle.  The  completion  of  an  unloading 
cycle  automatically  starts  the  next  unloading  cycle  until  the 
unloading  mission  of  the  B.U.F.  in  question  has  been  completed. 


* ' 

The  unloading  mission  of  the  Ath  B.U.F.  is  assumed  to  be  com¬ 
pleted  when  the  Ath  B.U.F.  is  not  needed  to  unload  any  mere 
payload  units. 


v)  Mission  of  the  ith  (i  =  1,2... I)  T.V. 


Step  1*  The  ith  T.V.  starts  from  its  base  having  as  its 
final  destination  W.A.I,  where  it  joins  queue  I. 


Step  2**  The  ith  T.V.  waits  in  queue  I  after  it  arrives  in 
W.A.I  until  the  appropriate  ship  unloading  area  is 
freed.***  (As  mentioned  earlier  at  the  start  of  the 
mission,  as  soon  as  the  Mother  Ship  is  moored,  all 
ship  unloading  areas  become  free.  Subsequently,  a 
ship  unloading  area  becomes  free  as  soon  as  the  T.V. 
that  is  being  served  alongside  the  ship  unloading 
area  in  question  is  unhooked  and  made  ready  to  com¬ 
mence  its  journey  to  W.A.II.) 


Step  3****  The  ith  T.V.  departs  from  W.A.I  for  the  Mother  Ship 
where  it  hooks  up  to  the  appropriate  ship  unloading 
area  as  soon  as  it  is  reached,  and  then  the  ith  T.V. 
is  made  ready  for  the  loading  operation.  The  T.V.  in 
question  departs  from  W.A.I  only  after  the  appro¬ 
priate  ship  unloading  area  is  free  to  receive  it. 


*Steps  1-13  will  not  be  executed  if  the  T.V.  in  question  breaks 
down  while  on  route  to  W.A.I  from  its  base.  The  ith  T.V.  is 
then  assumed  lost  from  our  system. 

**Steps  2-13  will  not  be  executed  if  the  T.V.  in  question 
breaks  down  while  waiting  in  W.A.I.  The  tth  T.V.  is  then  as¬ 
sumed  lost  from  our  system. 

***In  the  sequential  loading  mode  when  one  ship  unloading  area 
is  not  free,  then  all  ship  unloading  areas  are  considered  busy. 

****Steps  3-13  will  not  be  executed  if  the  T.V.  in  question 
breaks  down  while  on  route  from  W.A.I  to  the  appropriate  ship 
unloading  area.  The  ith  T.V.  is  then  assumed  lost  from  our 
system. 


Step  4  The  ith  T.V.  is  refuelled  (if  refueling  is  necessary) . 


Step  5  The  appropriate  payload  unit  is  secured  onto  the  ith 
T.V.  after  the  payload  unit  in  question  has  been  un¬ 
loaded  into  the  ith  T.V.  and  freed  from  the  appro¬ 
priate  S.U.F.,  and  after  the  S.U.F.  in  question  has 
been  made  ready  to  travel.  If  the  time  taken  in  this 
trip  for  the  ith  T.V.  to  secure  the  payload  unit  pre¬ 
viously  unloaded  (if  there  was  one)  onto  the  ith  T.V. 
is  less  than  the  time  taken  to  execute  Steps  2-5 
given  in  the  mission  of  the  kth  (1,2... K)  S.U.F. , 
then  the  ith  T.V.  waits  until  the  above  mentioned 
Steps  2-5  are  completed. 


Step  6  The  ith  T.V.  is  unhooked  from  the  ship  unloading  area 
in  question  and  is  made  ready  to  travel  after  the 
last  payload  unit  to  be  unloaded  into  the  ith  T.V.  in 
this  trip  of  the  T.V.  ir  question  has  been  properly 
secured.  (As  mentioned  earlier,  upon  completion  of 
this  step,  the  ship  unloading  area  in  question  be¬ 
comes  free.) 


Step  7*  The  ith  T.V.  starts  from  the  ship  unloading  area  in 

question  having  as  its  final  destination  W.A. II,  where 
it  joins  queue  II  after  it  has  been  made  ready  to 


_  travel. _ _ 

•Steps  7-13  will  not  be  executed  if  the  ith  T.V.  breaks  down 
while  on  route  from  the  ship  unloading  area  in  question  to  W.A. 
II.  The  payload  that  was  being  transported  on  this  trip  of  the 
ith  T.V.  and  the  T.V.  itself  are  then  assumed  lost  from  our 
system. 


I 


Step  8*  The  ith  T.V.  waits  in  queue  XI  after  it  arrives  at 
W.A.II,  until  the  appropriate  beach  unloading  area 
is  freed**.  (As  mentioned  earlier  at  the  start  of 
the  mission#  as  soon  as  the  Jlth  (l  =  1,2...L)  B.TJ.F. 
arrives  at  point  A,  the  Jlth  beach  unloading  area  be¬ 
comes  free.  Subsequently,  a  beach  unloading  area 
becomes  free  as  soon  as  the  T.V.  that  is  being 
served  at  the  beach  unloading  area  in  question  is 
made  ready  to  travel  again  for  W.A.I  or  its  base.) 

Step  9***  The  ith  T.V.  departs  from  W.A.II  for  the  shore, 

where  it  beaches  at  the  appropriate  beach  unloading 
area  as  soon  as  it  is  reached  and  then  the  ith  T.V. 
is  made  ready  for  the  unloading  operation.  The  T.V. 
in  question  departs  from  W.A.II  only  after  the  appro¬ 
priate  beach  unloading  area  is  free  to  receive  it. 


♦Steps  8-13  will  not  be  executed  if  the  T.V.  in  question 
breaks  down  while  waiting  in  W.A.II.  The  payload  that  was  be¬ 
ing  transported  on  this  trip  of  the  ith  T.V.  and  the  T.V.  it¬ 
self  are  then  assumed  lost  from  our  system. 

**In  the  sequential  unloading  mode,  when  one  beach  unloading 
area  is  not  free,  then  all  beach  unloading  areas  are  con¬ 
sidered  busy. 

***Steps  9-13  will  not  be  executed  if  the  T.V.  in  question 
breaks  down  while  on  route  from  W.A.II  to  the  appropriate 
beach  unloading  area.  The  payload  that  was  being  trans¬ 
ported  on  this  trip  of  the  ith  T.V.  and  the  T.V.  itself  are 
then  assumed  lost  from  our  system. 


Step  10  The  appropriate  payload  unit  is  released  from  the  ith 
T.V.  by  the  appropriate  B.U.F.  after  the  B.U.F.  in 
question  has  reached  the  tth  T.V.  Next,  the  appro¬ 
priate  payload  unit  is  secured  onto  the  B.U.F.  in 
question.  While  the  above  mentioned  B.U.F.  is  execut¬ 
ing  first  the  sixth  and  then  the  third  step  described 
in  the  mission  of  the  ith  (SL  -  lt2...L)  B.U.F.,  the 
ith  T.V.  waits  in  the  beach  unloading  area  in  ques¬ 
tion  until  the  above  mentioned  two  steps  are  com¬ 
pleted  . 

Step  11  The  £th  T.V.  is  made  ready  to  travel  after  the  last 
payload  unit  carried  by  the  ith  T.V.  in  this  trip 
of  the  T.V.  in  question  has  been  secured  onto  the 
appropriate  B.U.F.  (As  mentioned  earlier  upon 
completion  of  this  step  the  beach  unloading  area 
in  question  becomes  free.) 


Step  12* 


The  ith  T.V.  starts  from  the  beach  unloading  area 
in  question  having  as  its  final  destination  W.A.I, 
where  it  joins  queue  I  after  it  has  been  made  ready 
to  travel.  This  operation  is  performed  only  if  the 
transporting  mission**  of  the  ith  T.V.  has  not  been 
completed.*** 


Step  13****The  ith  T.V.  starts  from  the  beach  unloading  area 
in  question  after  it  has  been  made  ready  to  travel, 
having  as  its  final  destination  the  T.V. 's  base. 


Upon  completion  of  Step  13  given  above,  the  mission  of 
the  ith  T.V.  is  completed.  Please  note  that  when  Steps  2-12 
given  above  are  executed,  the  ith  T.V.  is  said  to  have  execut¬ 
ed  one  complete  trip.  The  completion  of  a  trip  automatically 
starts  the  next  trip  until  the  transporting  mission  of  the 

•  i 

T.V.  in  question  is  completed. 


♦Steps  12  and  13  will  not  be  executed  if  the  ith  T.V.  breaks 
down  while  on  route  from  the  beach  Unloading  area  in  question 
to  W.A.I.  The  ith  T.V.  is  then  assumed  lost  from  our  system. 

**The  transporting  mission  of  the  ith  T.V.  is  assumed  to  be 
completed  when  the  ith  T.V.  is  not  needed  to  transport  any 
more  payload  units. 

***In  our  mathematical  model  Step  12  is  always  executed  as 
long  as  there  is  even  one  payload  unit  on  board  the  Mother 
Ship.  However,  this  action  is  forfeited  for  all  the  T.V.  that 
happen  to  be  in  W.A.I  immediately  after  the  appropriate  T.V. 
departed  W.A.I  to  load  the  last  payload  unit(s)  from  the 
Mother  Ship. 

****Step  13  will  not  be  executed  if  the  T.V.  breaks  down  while 
on  route  from  the  beach  unloading  area  in  question  to  the 
T.V. 's  base.  The  ith  T.V.  is  then  assumed  lost  from  our 
system. 
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5.  Solution  Method 

For  the  reasons  given  in  Section  2 ,  the  method  of 
digital  simulation  is  selected  to  solve  the  problem  under  in¬ 
vestigation.  Because  of  the  size  of  the  problem,  the  use  of 
digital  computers  is  inevitable.  The  computer  program  for 
this  method  became  most  efficient  when  it  was  so  structured 
that  the  main  program  served  as  a  scheduler,  passing  control 
to  the  appropriate  subroutine,  which  simulated  the  appropriate 
event  at  the  time  when  control  was  passed.  Collectively,  the 
events  simulate  the  entire  mathematical  model  of  our  system. 

The  definition  of  the  events  and  of  the  scheduling  mechanism 
used  in  our  study  is  given  below. 

A)  Event  Definition 

The  events  selected  for  our  study  are: 

1.  i)  The  arrival  and  mooring  of  Mother  Ship, 

ii)  The  preparation  of  all  S.U.F.  for  the  unloading  opera¬ 
tion  and  their  arrival  at  the  appropriate  ship  unload¬ 
ing  areas, 

iii)  The  arrival  of  all  B.U.F.  at  point  A  on  the  beach  and 
their  preparation  for  the  unloading  operation,  and 

iv)  The  arrival  of  all  T.V.  in  W.A.I*. 

2.  a)  The  selection  of  a  T.V.  from  W.A.I,  or 
b)  The  selection  of  a  T.V.  from  W.A.II. 

*If  there  are  any  breakdown  considerations  involved  in  our  in¬ 
vestigation  for  this  part  of  the  trip  of  the  T.V.  in  question, 
then  they  must  be  taken  into  account. 


3.  a)  The  selection  of  a  S.U.F.  or 
b)  The  selection  of  a  B.U.F. 

\ 

4.  a)  The  loading  operation  alongside  the  Mother  Ship 

comprising 

i)  Steps  2-4  of  the  mission  of  the  ith  (i=l,2. .  .1)  T.V.* 

ii)  Steps  2-6  of  the  mission  of  the  fcth  (k-1, 2. . . K) 

S. U.F.**. 

iii)  Steps  5-7  of  the  mission  of  the  ith  (i=l,2. .. I)  T.V.* 
or 

b)  The  unloading  operation  at  the  beach  comprising 

i)  Steps  8  and  9  of  the  mission  of  the  ith  (i«*i,  2.  ..I) 

T. V.  * 

ii)  Steps  3-6. of  the  mission  of  the  ith  (1=1,2. . .L)  B.U.F. 

iii)  Steps  10-12  of  the  mission  of  the  ith  (i=l,2. . .1) 

T.V.  * 

5.  The  closing  of  the  entire  operation  ***  by  executing 

i)  Event  2b, 
ii)  Event  3b, 

# 

If  there  are  any  breakdown  considerations  involved  in  our  in¬ 
vestigation  for  this  (these)  part(s)  of  the  trip  of  the  T.V. 
in  question,  then  they  roust  be  taken  irto  account. 

**Please  note  that  as  soon  as  the  ftth  (k=l, 2. . . K)  S.U.F.  un¬ 
loads  its  entire  payload,  step  7  of  the  mission  of  the  kth 
(k=l,2...K)  S.U.F.  is  executed  for  the  S.U.F.  in  question. 

***This  event  is  executed  only  when  there  are  no  more  payload 
units  to  be  unloaded  from  the  Mother  Ship. 


iii)  Steps  8  and  9  of  the  mission  of  the  ith  (i=l ,  2. . . I) 
T.V.*  for  all  T.V.  in  W.A.II, 

iv)  Steps  3-6  of  the  mission  of  the  Ith  ( 1=1, 2. . .  L) 

B.U.F. , 

v)  Steps  10  and  11  of  the  mission  of  the  ith  (i=l,2...I) 
T.V.*  for  all  T.V.  in  W.A.II/ 

vi)  Step  3  of  the  mission  of  the  Mother  Ship, 

vii)  Steps  7  and  8  of  the  mission  of  the  Jtth  (1=1,2. .  .L) 
B.U.F. ,  and 

viii)  Step  13**  of  the  mission  of  the  ith  (i=l,2. . .1)  T.V.* 


B)  Scheduling  Mechanism  Definition 

The  considerations  leading  to  the  definition  of  the 
scheduling  mechanism  that  is  to  be  used  by  the  main  program 
for  scheduling  the  appropriate  event  at  any  given  time  are 
given  below. 

At  the  outset  of  our  investigation  it  was  found  essential 
that  our  model  should  be  provided  with  a  simulation  clock, 
which  is  to  be  used  to  record  the  start  of  an  event  execution, 
as  this  information  was  considered  to  be  useful  output.  Next, 


*If  there  are  any  breakdown  considerations  involved  in  our  in¬ 
vestigation  for  this  part  of  the  trip  of  the  T.V.  in  question, 
then  they  must  be  taken  into  account. 

**The  action  taken  by  executing  step  12  of  the  mission  of  the 
ith  (i=l, 2. . . I)  T.V.  is  forfeited  for  all  the  T.V.  that  are 
in  W.A.I  when  the  execution  of  the  fifth  event  commences. 


* 


it  was  noted  that  as  our  simulation  is  event  structured, 
the  simulation  clock  could  also  be  used  as  the  scheduling 
mechanism  if 

i)  at  the  start  of  the  entire  operation  the  clock  was 
initialized  to  a  zero  reference  time, 

ii)  at  subsequent  times  the  clock  was  always  updated  to 
show  the  starting  time  of  the  last  event  executed,  and 

iii)  the  time  at  which  an  event  is  terminated  is  continuously 
updated . 

This  is  so  because  we  may  now  schedule  the  next  event  correctly 
by  simply  selecting  the  earliest  available  event  after  that 
just  executed,  as  recorded  by  the  simulation  clock.  The  above 
description  defines  the  scheduling  mechanism  employed  in  this 
study  for  scheduling  the  event  that  is  to  be  executed  at  any 
given  time. 

In  our  mathematical  model,  once  the  correct  event  is 
decided  by  using  the  above  mentioned  method,  the  action  taken 
is  as  follows: 

Control  is  passed  to  the  appropriate  subroutine,  which 
executes  the  event  in  question,  updates  all  the  relevant 
variables  and  the  simulation  clock,  and  then  returns  con¬ 
trol  to  the  scheduler,  which  repeats  the  above  procedure 
until  the  fifth  event  is  executed,  terminating  the 
analysis  of  the  run  under  investigation. 


From  the  above,  it  follows  that  because  the  simulation 
clock  advances  only  in  discrete  jumps  (including  the  zero 
jump) ,  then  we  may  collect  any  information  relevant  to  the 
problem  by  examining  the  simulation  only  at  the  discrete 
times  recorded  by  the  simulation  clock.  This  is  so  because 
no  event  may  start  at  any  other  time,  and  as  the  entire  simu¬ 
lation  is  represented  by  events,  no  additional  information 
may  be  obtained  by  examining  the  simulation  more  frequently. 

The  above  discussion  completes  our  description  of  the 
event  structure  of  our  mathematical  model  and  the  scheduling 
mechanism,  the  two  most  important  aspects  of  our  method.  All 
other  aspects  of  our  method,  such  as  random  number  generation, 
tests  for  random  number  generators,  etc.,  are  not  covered 
here  as  their  treatment  may  be  found  elsewhere. 

To  permit  the  reader  to  fully  understand  the  simulation 
process,  a  brief  description  of  the  computer  program  is  given 
below,  together  with  a  general  flow  chart*,  which  will  also 
serve  to  show  the  reader  how  to  introduce  new  use  strategies 
when  this  is  found  desirable. 


For  the  sake  of  simplicity  of  presentation,  the  general  flow 
chart  does  "ot  include  the  breakdown  considerations. 


0 


92 


C.  Computer  Program  Description. 

MAIN 

The  MAIN  program  in  our  simulation  serves  as  a  scheduler,  and 

operates  as  follows: 

1.  Reads  the  input  specifying  the  number  of  cases  to  be 
processed  by  this  computer  run. 

2.  Passes  control  to  Subroutine  INPUT  to  read  the  input. 

3.  Passes  control  to  Subroutine  INOUT  to  print  the  input. 

4.  Initializes  the  model. 

5.  Passes  control  to  Subroutine  BEGIN  to  simulate  the  arrival 
of  our  resources  in  the  theatre  of  operations. 

6.  Passes  control  to  Subroutine  ASLTVA  or  ASLTVB  which  in  turn 
passes  control  to  Subroutine  BSLTVA  or  BSLTVB . 

7.  Passes  control  to  Subroutine  SLSCA  or  SLSCB  or  to  Sub¬ 
routine  SLBCA  or  SLBCB,  depending  on  the  outcome  of  the 
above  step. 

8.  Passes  control  to  Subroutine  LOAD  or  UNLOAD  or  repeats 
steps  6  and  7  given  above,  depending  on  the  outcome  of  the 
two  steps  just  mentioned. 

9.  Passes  control  to  Subroutine  FIN,  which  in  turn  passes  con¬ 
trol  to  Subroutines  BSLTVA  or  BSLTVB  and  Subroutines  SLBCA 
or  SLBCB,  until  all  T.V.  in  W.A.II  are  processed.  Sub¬ 
routine  FIN  then  simulates  the  departure  of  our  resources 
from  the  theatre  of  operations  at  the  conclusion  of  our 
mission,  and  finally, 


9  3 
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10.  If  no  more  runs  are  to  be  executed  for  the  case  under 

investigation,  control  is  passed  to  Subroutine  STATIC  to 
compute  the  necessary  statistics.  If  more  runs  are  to  be 
executed,  steps  4-9  given  above  are  repeated  until  the 
number  of  runs  executed  is  equal  to  that  prespecified. 

It  should  be  noted  that  the  program  is  so  arranged  as  to 
permit  the  processing  of  as  many  cases  as  is  desired  by  simply 
providing  the  necessary  input  information.  In  addition,  it 
should  be  noted  that  during  the  execution  of  any  run,  as  soon 
as  all  transfer  vehicles  are  found  to  be  malfunctioning,  the 
processing  of  this  run  is  terminated  and  the  processing  of  the 
next  one  (if  any)  is  automatically  started. 

INPUT 

This  Subroutine  reads  all  the  input  pertaining  to  the  case 
under  investigation. 

The  input  for  each  case  comprises: 

1.  The  specification  of  the  number  of  runs. 

2.  The  specification  of  the  units  of  the  variables. 

3.  The  specification  cf  the  variables  defining  the  payload 
characteristics  and  its  allocation. 

4.  The  specification  of  the  variables  defining  the  Mother  Ship 
characteristics . 

5.  The  specification  of  the  variables  defining  the  S.U.F.,  the 


B.U.F.  and  the  T.V. 


In  addition/  the  input  specifies  a)  the  loading  mode  at  the 
Mother  Ship#  b)  the  unloading  mode  at  the  beach,  and  c)  the 
use  strategies  to  be  used  at  each  of  the  decision  nodes. 
Furthermore,  the  input  specifies  whether  a  process  is  deter¬ 
ministic  or  stochastic.  If  the  process  is  stochastic,  the 
input  designates  the  type  of  random  number  generator  to  be 
used,  its  seed,  and  the  range  of  variation  of  the  random  num¬ 
bers.  Finally,  the  input  indicates  whether  malfunctioning 
considerations  are  to  be  included  in  our  analysis  and,  if  so, 
supplies  the  necessary  information. 

Please  note  that,  as  was  already  mentioned  in  Section  3,  the 
above  input  will  be  appropriately  augmented  in  order  to  permit 
the  specification  of  the  statistical  considerations  and  of  the 
antithetic  variance  concept.  These  considerations  and  this 
concept  will  be  developed  in  a  separate  report,  which  will 
follow  the  present  one. 


INPUT 

This  Subroutine,  when  developed,  will  print  out  all  the  input 
pertaining  to  the  case  under  investigation.  This  subroutine 
has  not  yet  been  developed,  as  it  is  not  needed  for  the  comput¬ 
er  program  testing.  Once  the  computer  program  is  tested  and 
its  final  form  is  decided,  then  the  INOUT  subroutine  will  be 
written.  In  this  way,  any  changes  and  additions  in  the  present 
input  will  be  directly  incorporated  into  INOUT  with  minimum 


labor. 


RANDU 


This  Subroutine  computes  uniformly  distributed  random  real 
numbers  between  zero  and  one.  It  is  of  the  congruential  type 
and  is  modeled  after  the  subroutine  listed  in  IBM's  Scientific 
Subroutine  Package*.  With  a  seed  equal  to  65539,  it  has  a 
cycle  length  of  229  terms  and  it  satisfies  all  the  usual  tests 
for  randomness. 

BEGIN 

This  Subroutine  simulates  the  first  event,  described  at  the 
beginning  of  this  section.  If  any  breakdown  considerations  are 
to  be  included  in  this  investigation  for  this  part  of  the  sim¬ 
ulation,  all  T.V.  found  to  be  malfunctioning  are  refused  entry 
into  W.A.I. 

ASLTVA 

This  Subroutine  defines  the  first  of  the  strategies  incorpor¬ 
ated  in  the  model  for  selecting  T.V.  from  W.A.I**.  It  is 
basically  a  first  come,  first  served  strategy  and  is  oriented 
towards  volume-limited  transfer  vehicles.  The  rules  utilized 
by  Subroutine  ASLTVA  for  the  above  mentioned  T.V.  selection 
have  already  been  given  in  Section  4  (see  T.V.  use  strategy  A 
for  W.A.I) . 

^Seethe  bibliography  at  the  end  of  this  report. 

♦♦Compare  with  event  2a,  described  at  the  beginning  of  this 
section. 


ASLTVB 


This  Subroutine  defines  the  second  of  the  strategies  incor¬ 
porated  in  the  model  for  selecting  T.V.  from  W.A.I*.  It  is 
again  basically  a  first  come,  first  served  strategy  but  is 
oriented  towards  weight-limited  transfer  vehicles.  The  rules 
utilized  by  Subroutine  ASLTVB  for  the  above  mentioned  selec¬ 
tion  have  already  been  given  in  Section  4  (see  T.V.  use 
strategy  B  for  W.A.I). 


BSLTVA 

This  Subroutine  defines  the  first  of  the  strategies  incorpor¬ 
ated  in  the  model  for  selecting  T.V.  from  W.A.II**.  It  is 
basically  a  first  come,  first  served  strategy  and  is  oriented 
towards  volume-limited  transfer  vehicles.  The  rules  utilized 
by  Subroutine  BSLTVA  for  the  above  mentioned  T.V.  selection 
have  already  been  given  in  Section  4  (see  T.V.  use  strategy  A 
for  W.A.II) . 

BSLTVB 

This  Subroutine  defines  the  second  of  the  strategies  incorpor¬ 
ated  in  the  model  for  selecting  T.V.  from  W.A.II?*  It  is  again 
basically  a  first  come,  first  served  strategy  but  is  oriented 
towards  weight-limited  transfer  vehicles.  The  rules  utilized 

Compare  with  event  2a,  described  at  the  beginning  of  this 
section. 

**Compare  with  event  2b,  described  at  the  beginning  of  this 
section. 


by  Subroutine  BSLTVB  for  the  above  mentioned  T.V.  selection 

\ 

have  already  been  given  in  Section  4  (see  T.V.  use  strategy  B 
for  W. A. II). 

SLSCA 

This  Subroutine  defines  the  first  of  the  strategies  incorpor¬ 
ated  in  the  model  for  selecting  S.U.F.*  It  is  a  first  come, 

first  served  strategy.  The  rules  utilized  by  Subroutine  SLSCA 

\ 

for  the  above  mentioned  S.U.F.  selection  have  already  been 
given  in  Section  4  (see  S.U.F.  use  strategy  A). 

SLSCB 

This  Subroutine  defines  the  second  of  the  strategies  incorpor¬ 
ated  in  the  model  for  selecting  S.U.F.*.  It  is  basically  a 
first  come,  first  served  strategy.  The  rules  utilized  by  Sub¬ 
routine  SLSCB  for  the  above  mentioned  S.U.F.  selection  have 
already  been  given  in  Section  4  (see  S.U.F.  use  strategy  B) . 


LOAD 

This  Subroutine  simulates  event  4a,  described  at  the  beginning 

of  this  section.  If  any  breakdown  considerations  are- to  be 

included  in  this  investigation  for  this  part  of  the  simulation, 

all  T.V.  found  to  be  malfunctioning  are  refused  departure  from 

W.A.I,  arrival  at  the  ship  unloading  area  in  question,  or  entry 

into  W.A.II,  depending  upon  where  the  breakdown  occurred. 

•Compare  with  event  3a,  described  at  the  beginning  of  this  sec¬ 
tion. 


SLBCA 


This  Subroutine  defines  the  first  of  the  strategies  incorpor¬ 
ated  in  the  model  for  selection  B.U.F.*  It  is  a  first  come, 
first  served  strategy.  The  rules  utili.zed  by  Subroutine  SLBCA 
for  the  above . mentioned  B.U.F.  selection  have  already  been 
given  in  Section  4  (see  B.U.F.  use  strategy  A). 


SLBCB 

This  Subroutine  defines  the  second  of  the  strategies  incorpor¬ 
ated  in  the  model  for  selecting  B.U.F.*  It  is  basically  a 
first  come,  first  served  strategy.  The  rules  utilized  by  Sub¬ 
routine  SLBCB  for  the  above  mentioned  B.U.F.  selection  have 
already  been  given  in  Section  4  (see  B.U.F.  use  strategy  B) . 

UNLOAD 

This  Subroutine  simulates  event  4b,  described  at  the  beginning 
of  this  section.  If  any  breakdown  considerations  are  to  be 
included  in  this  investigation  for  this  part  of  the  simulation 
all  T.V.  found  malfunctioning  are  refused  departure  from 
W.A.II,  arrival  at  the  beach  unloading  area  in  question,  or 
entry  into  W.A.I,  depending  upon  where  the  breakdown  occurred. 


♦Compare  with  event  3b,  described  at  the  beginning  of  this  sec 
£ion. 
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FIN 

This  Subroutine  simulates  the  fifth  event,  described  at  the 
beginning  of  this  section.  If  any  breakdown  considerations 
are  to  be  included  in  this  investigation  for  this  part  of  the 
simulation,  all  T.V.  found  malfunctioning  are  refused  depart¬ 
ure  from  W.A.II,  arrival  at  the  beach  unloading  area  in  ques¬ 
tion,  or  arrival  at  their  bases,  depending  upon  where  the 
breakdown  occurred. 

STATIC 

Provisions  have  been  made  to  allow  the  introduction  of  sub¬ 
routines  which,  when  developed,  will  gather  the  sample  points 
generated  by  our  analysis,  compute  the  statistics  of  interest, 
and  output  all  quantities  necessary  for  the  final  evaluation 
of  alternative  strategies.  The  development  of  these  subrou¬ 
tines  is  deferred  until  the  form  of  the  desired  output  has 
been  finalized,  so  that  any  changes  and  additions  from  the 
present  output  will  be  directly  incorporated  into  STATIC  with 
minimum  labor. 

j1 

The  above,  together  with  Fig.  5-1  (see  next  page)  completes 
our  discussion  of  the  solution  method.  We  now  proceed  with 
the  discussion  of  the  problem  solution. 


SLBCA 


SLBCB 


Introduce  New 
.Strategies  Here. 


Fig.  5-1  (continued) 
General  Flew  Chart 


END 


103 


6.  Problem  Solution 

The  data  setup  for  the  computer  program  is  shown  in  Figs. 
6.1  -  6.17*,  and  must  be  followed  if  the  computer  program  is 
to  give  a  meaningful  solution  to  the  problem. 

Because  the  development  of  the  concept  of  antithetic  var¬ 
iance  is  deferred  until  the  next  report,  it  is  more  profitable 
for  us  to  delay  the  preparation  of  the  input  and  the  solution 
of  the  problem  until  then.  It  therefore  follows  that  the 
last  topic  of  the  solution  process,  namely,  the  evaluation  of 
results,  must  also  be  deferred  until  then. 

The  above  completes  the  main  body  of  this  report.  In 
the  next  and  final  section,  future  recommendations  drawn  from 
our  analysis  are  given.  Finally,  in  Appendix  A  of  this 
report,  the  computer  program**  is  listed  for  easy  reference. 


*The  variables  controlling  the  statistical  considerations  and 
the  antithetic  variance  concept  will  be  introduced  in  the 
next  report. 

**Please  note  that  the  computer  program  has  been  debugged,  but 
its  logic  has  not  been  exhaustively  checked  because  it  is 
more  economical  to  await  the  development  of  the  concept  of 
the  antithetic  variance. 


*1.  General  Data  Setu 


Input  Data  Describing  tne  Breakdown  Consider¬ 
ations  for  the  j^th  Case  (for  a  more 
detailed  breakdown  refer  to  Fig.  6-17) 


Input  Data  Describing  the  Stochastic 
Behavior  of  the  j^th  Case  (for  a  more 
detailed  breakdown  refer  to  Fig.  6-13) 

Input  Data  Indicating  the  Stochastic 

Behavior  of  the  j^th  Case  (for  a  _ _ 

more  detailed  breakdown  refer  to 
Fig.  6-9) 


WGHT  VOL  TIME 


ISUD  IBUD 
IWA1SL  IWA2SL 
ISCSL  IBCSL 


Input  Data  for  Time  Character^ 
istics  (for  a  more  detailed 
breakdown  refer  to  Fig.  6-5) 


Input  Data  for  T.V.'s  Payload 
Characteristics  (for  a  f 

more  detailed  breakdown  / 
refer  to  Fig.  6-4) 


I » K ,  L 


7 


-X 


v 


y 


y 


Input  Data  for 
Payload  De¬ 
scription  (for 
a  more  detailed 
breakdown  refer 
to  Fig.  6-3) 


Number  of  T.V.,  S.U.F. 
and  B.U.F.  involved  in 
the  j^th  case 

.Number  of  runs  to  be  processed 
for  the  j-^th  case 


Fig.  6-2.  Input  Data  Setup  for  the  j^th  Case  ( j^=l,2. . .NCASES) 


l| 
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Fig.  6-4 

Input  Data  Setup  for  the  T.V.'s  Payload  Characteristics 


Input  Data  Setup  for  the  B.U.F.'s  Time  Characteristics 


Input  Data  Indicating  the  Stochastic 
^Behavior  of  the  T.V.  (for  a  more 
detailed  breakdown  refer  to  Fig.  6-12) 


Input  Data  Indi¬ 
cating  the  Sto¬ 
chastic  Behav-  / 
ior  of  the  / 

B.U.F.  - w4— 

(for  a  more  _ 

detailed  / 

breakdown  / _ 

refer  to 

Fig-  6-U|  '  A  T 


IDMA, ID1 , ID2 


Input  Data  Indicating  the-- 
Stochastic  Behavior  of  the 
Mother  Ship 


'Input  Data" Indicat 
ing  the  Stochastic 
Behavior  of  the 
S.U.F.  (for  a  tore 
detailed  breakdown 
refer  to  Fig.  6-10 


Input  Data  Setup  for  the  Indication  of  the 


Stochastic  Behavior  of  the  j'^th  Case 
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IDSC9,1  IDSC9,2 


IDSC8,1  IDSC8,2 


IDSC7,1  IDS  C  7 , 2 


IDSC_  .  IDSC.  . 
5,1  5,2 


IDSC4,1  IDSC4,2 


me Tivsc- . 

J  ,  i.  J  ,  z 


IDSC2,1  IDSC2,2 


IDSC1,1  IDSC1,2 


IDSC 


9.K 


IDSC 


8 ,  K 


IDSC 


7  >  K 


IDSC 


5  ,Kj 


IDSC 


4,K 


IDSC3,K 


IDSC 


2  ,  K 


IDSC 


1,K 


IDBC8,1  IDBC8,2 

[  IDBC7,1  IOBC7,2  * 

IDBC6,1  IDBC7~ 

IDBC4,1  IDBC4,2  . 

IDBC3,1  IDBC3,2  . 

idbc2,i  idbc7; .......... 

IDBC1,1  IDBC1,2  . 


IDTV 


Input  Data  Setup  for  the 
Indication  of  the  Stochastic 
Behavior  of  the  T.V. 


Input  Data  Describing 
the  Stochastic  Behavior 
of  the  B.U.F.  (for  a  more 
detailed  description 


Input  Data  Describing 
the  Stochastic  Behavior 
of  the  T.V.  (for  a  more 
detailed  description 


If  ID2-1,  remove  this  card  from  the  input 


Fig.  6-13 

Input  Data  Setup  for  the  Description  of 
the  Stochastic  Behavior  of  the  jjth  Case 


INSC9  k  AMINSCg  K 

/  amaxsc9>k 


INSC9  2 AMINSC9  2 
’  AMAXSCj’j 

I N  S  C  q  i  AMINSC,  , 

'  AMAXSC9  j 

INSCg  K  AMINSCg  K 

/  AMAXSC.  „ 
/  8  ,  K 


INSCg  ,  AMINSCg  , 

’  AMAXSC * ’ x 

5C7  ,K  AMINSd71c 

/  *maxsc7k 


I  NSC  7  ,  AMINSC7  , 
’  AMAXSC 7*2 

fSC  .  ,  AMINSC, 

’  AMAXSC  7 ’ , 


/  AMAXSC,  v 
/  5  ,  X 


'3  2  AMI.Otj  2 

’  amaxsc5’2 
l  aminsc51 

/  AMAXSC c 


If  IDSCj2>k  (j2- 

1,2.  .  .5,7,8,9,10; 
k“1 , 2 . . . K)  equals  1 , 
remove  the  (jj.lOth 
card  from  the  input 


INSC 7  K  AMINSCj  K 

•  AMAXSC. 


INSC,  ,  AAlNSC,  , 
AMAXSCj  2 
NSC,  ,  AMINSC,  , 

'  AMAXSC  2  *  ^ 

"l.K  AMINSC1,X  ’ 

/  AMAXSC,  „ 


INSC12 

AMINSC,  j 
AMAXSC} | 2 

INSC 1 , 1 

AMINSCg  x 

AMAXSC 1  : 

Fig.  6-14 


Input  Data  Setup  for  the  Description  of 
the  Stochastic  Behavior  of  the  S.U.F. 


* 
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7 .  Future  Recommendations 

Before  concluding  this  report,  I  would  like  to  take  the 
opportunity  of  discussing  what  I  envision  as  Phase  II  (the 
next  report)  of  this  project. 

First,  the  development  of  the  antithetic  variance  concept 
for  general  congestion  models.  This  not  only  applies  to  our 
case  but  also  to  any  general  congestion  problem.  The  time¬ 
saving  element  will  be  the  major  reward  of  such  an  analysis, 
as  it  is  anticipated  that  the  simulation  will  be  made  twice  as 
efficient.  (The  comparison  is  made  between  a  simulation 
utilizing  independent  samples  and  one  utilizing  negatively 
correlated  samples.)  Secondly,  the  application  of  the  concept 
of  antithetic  variance  to  our  model.  Thirdly,  the  development 
of  the  subroutine  INOUT  and  of  the  subroutine  package  STATIC, 
and  fourthly  and  finally,  recommendations  for  a  possible  third 
stage  which  will  include  the  graphic  representation  of  our 
simulation  as  it  is  executed.  This  I  believe  to  be  very  help¬ 
ful,  because  the  user  can  then  follow  the  model  visually  and 
impose  logical  changes  by  assuming  manual  control  whenever  it 
becomes  obvious  that  a  previously  selected  strategy  is  inappli¬ 
cable.  This  has  been  tried  in  the  Civil  Engineering  Depart¬ 
ment  of  M.I.T.  for  a  car  assignment  simulation,  with  relatively 
high  success. 
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